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Equation  (3.23)  (and  the  equation  following): 


tern  should  read; 
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The  right  hand 
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Equation  (3.Zh)  -  Replace  a  with  ^  so  that  the  expression 
should  read: 
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Equation  (B.3)  —  should  read  as  follows: 
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ABSTRACT 


>• 


This  work  is  concerned  with  tunable  light  scattering  from 
optical  modes  in  crystals.  The  motivation  for  this  study  is  the 
generation  of  tunable  optical  and  infrared  emission,  an  improved 
understanding  of  the  scattering  process,  and  a  study  of  material 
properties . 

The  process  involved  is  a  combination  of  the  Raman  and  para¬ 
metric  effects  in  a  crystal  with  infrared  absorbing  optical  modes. 

The  input  (pump)  photon  excites  a  lattice  vibration  and  a  signal 
photon  is  emitted  at  ths  difference  frequency  between  the  pump  photon 
and  the  lattice  vibration.  Simultaneously,  idler  radiation  is  gener¬ 
ated  at  the  vibration  frequency  due  to  the  infrared  activity  of  the 
scattering  optical  mode.  The  dispersion  of  the  idler  propagation 
characteristic  produces  tuning  of  the  signal  and  idler  radiation  as 
the  angle  between  the  pump  and  signal  beam  is  varied. 

A  classical  analysis  of  the  stimulated  process  is  presented. 

The  analysis  uses  an  energy  density  formulation  and  follows  the 
approach  of  Henry  and  Garrett.  The  contribution  of  the  present  work 
in  this  regard  is  to  extend  their  results  to  multimode  materials 
and  to  derive  a  gain  expression  which  is  valid  for  all  idler  fre¬ 
quencies.  A  quantum  mechanical  analysis  of  the  spontaneous  scatter¬ 
ing  process  is  given.  Expressions  for  the  scattered  signal  and  idler 
powers  are  derived  and  techniques  for  computing  the  material  nonlinear 
coefficients  are  presented.  Since  the  experiments  performed  in  this 
study  involved  interactions  with  the  transverse  optical  modes  of  lithium 
niobate  (LiNbO,.),  the  results  of  these  calculations  are  applied  in 
detail  to  this  case. 


The  most  significant  experimental  results  of  this  work  were 
obtained  in  stimulated  scattering  from  the  2^8  cm”1  mode  in  LiNbU^. 

A  1  MW  Q-swito'ned  ruby  laser  was  used  as  a  pump,  and  tunable  stimulated 
signal  and  idler  outputs  were  observed  without  the  ut  »  of  an  external 
resonator.  The  crystal  surfaces  were  polished  flat  and  parallel  to 
each  other  to  provide  for  several  passes  of  the  signal  radiation 
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through  the  pump  beam.  Tuning  was  obtained  simply  by  rotating  the 
normal  to  the  crystal  surfaces  with  respect  to  the  pump  beam.  Con¬ 
version  efficiencies  of  the  pump  to  the  signal  of  greater  than  5C$ 
were  measured.  While  the  tuning  range  of  the  signal  was  relatively 
small  (6965  8  -  70^5  ft),  the  idler  wavelength  was  tuned  from  50  |i 
to  250  n,  with  a  peak  power  output  of  about  5  watts.  The  experimental 
results  were  found  to  be  in  reasonable  agreement  with  predicted  values. 

Spontaneous  scattering  measurements  were  also  performed.  A 

significant  part  of  the  experimental  effort  in  this  area  was  devoted 

to  the  construction  of  an  argon  laser  and  the  assembly  of  a  phase- 

sensitive  detection  system.  This  scattering  test  set  is  adaptable 

to  a  wide  variety  of  spontaneous  scattering  experiments.  It  was 

used  to  investigate  the  properties  of  the  spontaneously  scattered 

signal  power  which  is  generated  in  scattering  from  the  248  cm  ^  mode 

in  LiNbO-,.  The  frequency  and  magnitude  of  the  signal  were  measured 

3 

as  a  function  of  the  scattering  angle.  Again,  the  results  were  found 
to  agree  closely  with  calculations. 
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CHAPTER  I 


INTRODUCTION 

A.  OBJECT  OF  THIS  WORK 

The  object  of  this  work  is  the  study  of  tunable  light 
scattering  from  optical  modes  in  crystals.  The  motivation  for 
this  effort  is  the  generation  of  tunable  optical  and  infrared 
emission,  an  improved  understanding  of  the  scattering  process,  and 
a  study  of  material  properties .  Spontaneous  and  stimulated  scat¬ 
tering  are  analyzed  theoretically  and  the  results  of  these  calcula¬ 
tions  are  compared  to  experimental  observations.  The  analysis  and 
experiments  are  specifically  concerned  with  interactions  with  the 
transverse  optical  modes  of  lithium  niobate  (LiNbOy, 

B.  HISTORICAL  BACKGROUND 


The  pioneering  study  of  this  type  of  light  scattering  was  done 

in  1928  by  Raman  and  Krishnan.1  They  observed  that  a  small  fraction 

of  the  light  transmitted  through  liquids  and  gases  emerged  with  a 

frequency  shifted  from  that  of  the  incident  source  by  the  vibrational 

resonant  frequency  of  the  scattering  medium.  This  phenomenon  is 

called  the  Raman  effect.  The  stimulated  Raman  effect  was  discovered 

2 

unexpectedly  in  1962  by  Woodbury  and  Ng,  in  the  output  of  a  Q- 


switched  ruby  laser.  This  discovery  was  quickly  followed  by  the 

5-6 

observation  of  stimulated  Raman  emission  from  other  materials  and 

5  7-10 

the  theoretical  explanation  of  the  effect  was  soon  developed,  ’ 


C.  TUNABLE  RAMAN  SCATTERING 


Raman  and  parametric  scattering  are  examples  of  processes  in 
which  an  input  wave  (pump),  passing  through  a  material,  generates 
waves  at  new  frequencies  (signal  and  idler).  Conservation  of  energy 


requires  that  oo^  =  +  un  ,  where  a>  is  the  frequency  and  p  , 

s  ,  and  •  refer  to  pump,  signal  and  idler  respectively.  Materials 
which  exhibit  this  property  are  in  a  class  of  materials  which  is 
termed  nonlinear. 

In  parametric  scattering,  the  idler  wave  is  purely  electro¬ 
magnetic  in  nature  and  the  nonlinearity  which  produces  the  signal  polari¬ 
zation  is  due  solely  to  the  motion  of  electrons  in  the  material,  since, 
in  general,  the  idler  frequency  is  far  above  the  highest  vibrational 
resonant  frequency.  This  process  produces  radiation  at  both  signal 
and  idler. 

In  nont unable  Raman  scattering,  the  idler  frequency  corresponds 
to  a  vibrational  resonance  of  the  nonlinear  material.  The  idler  wave 
consists  of  a  mechanical  vibration  and  the  nonlinearity  associated  with 
the  signal  polarization  is  due  to  ionic  motion.  Therefore,  radiation 
is  generated  at  only  the  signal  frequency.  Vibrations  which  exhibit 
this  effect  are  termed  Raman  active. 

If  separation  of  electrical  charge  occurs  during  the  optical 
vibration  of  a  crystal  lattice,  a  dipole  moment  at  the  vibration  fre¬ 
quency  is  generated.  The  material  can  then  both  absorb  and  generate 
electromagnetic  radiation  at  this  frequency  and  is  called  infrared 
active.  (  The  optical  vibrations  are  at  infrared  frequencies.)  The 
dispersion  curve  for  an  infrared  active  crystal  is  strongly  affected 
by  this  coupling  between  electromagnetic  radiation  and  the  optical  modes 
of  the  material.1'1"  The  dispersion  curve,  which  is  shown  in  Fig.  1.1, 

is  partly  phonon-like,  partly  photon-like,  and  partly  phonon-photon- 

12  13 

like.  The  quanta  associated  with  the  mixed  system  are  sometimes 

l4 

called  polantons. 

The  possibility  of  generating  tunable  radiation  from  materials 
which  are  simultaneously  infrared  and  Raman  active  was  suggested  inde¬ 
pendently  by  Loudon1^  and  Pantell  and  Puthoff.1^  Briefly  stated,  the 
tuning  is  accomplished  by  controlling  the  idler  propagation  constant 
k  ,  which,  in  turn,  determines  the  idler  frequency  through  the  disper¬ 
sion  curve.  As  the  propagation  constant  is  varied,  the  idler  frequency 
moves  from  the  phonon-like  region  of  the  dispersion  curve  (near  the 
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resonant  frequency)  to  a  mixed  region,  where  the  idler  wave  is  both 
mechanical  and  electronic gnetic  in  nature.  In  this  region,  both  elec¬ 
tronic  and  vibrational  nonlinearities  can  contribute  to  the  process, 
and  radiation  at  both  signal  and  idler  frequencies  is  generated.  The 
idler  propagation  constant  is  varied  by  adjusting  the  angle  between 
the  pump  and  signal  beams  which,  through  wave  vector  conservation, 
controls  the  idler  propagation  constant.  Thus,  tunable  signal  and 
idler  outputs  can  be  obtained  in  interactions  with  crystal  vibrations 
which  are  Raman  and  infrared  active . 

D.  SUMMARY  OF  RESEARCH  EFFORT 


At  the  inception  of  the  experimental  program,  stimulated 

scattering  from  LiNbO^  had  not  as  yet  been  observed.  Therefore,  the 

study  of  the  tunable  idler  power  generated  in  spontaneous  scattering 

from  the  628  cm  optical  mode  in  LiNbO,  was  selected  as  an  initial 

5  17 

goal.  This  was  a  continuation  of  the  work  of  Puthoff,  et  al.,  who 

had  examined  the  signal  power  generated  in  this  scattering  process. 

The  detection  of  the  idler  would  have  been  significant  for  two  reasons. 

First,  it  would  have  provided  a  simple  tunable  source  in  the  l6-20  " 

wavelength  region.  Second,  its  detection  would  have  verified  for  the 

1.5  18  21 

first  time,  the  theoretical  prediction  of  its  existence.  ’ 

For  the  purpose  of  the  experiment,  an  argon  laser  was  constructed 
for  use  as  a  pump,  and  a  phase- sensitive  detection  system,  incorpora¬ 
ting  a  Ge:Cu  infrared  detector,  was  assembled.  The  laser  and  detec¬ 
tion  system  form  a  test  set  for  spon:  neous  scattering  measurement 
which  is  useful  for  a  broad  range  of  such  measurements,  with  only 
minor  modifications. 

During  the  course  of  this  work,  however,  the  first  observation 

22 

of  stimulated  scattering  from  LiNbO^  was  reported  and  was  quickly 
followed  by  the  detection  of  tunable  stimulated  outputs  (using  an 
external  resonator)  from  the  same  material.  The  effort  on  this 


(2) 


program  was  therefore  shifted  to  an  investigation  of  the  stimulated 

Raman  effect.  In  the  course  of  these  experiments,  tunable,  stimulated, 

optical  scattering  from  the  248  cm  mode,  which  did  not  require  an 

external  resonator,  was  observed.  In  addition,  the  idler  radiation 

24 

accompanying  this  process  was  detected. 

Finally,  measurements  of  the  spontaneous  scattering  from  the 
248  cm  mode  were  performed. 

Theoretical  calculations  were  carried  out  in  conjunction  with 
the  experimental  effort.  The  results  of  these  calculations  were  in 
good  agreement  with  the  experimental  results  for  both  the  stimulated 
and  spontaneous  cases. 

In  summary,  the  highlights  of  this  work  are: 

(l)  The  elaboration  and  extension  of  the  theory  of  stimulated 
and  spontaneous  scattering  and  its  application  to  the 
A^- symmetry  modes  of  LiNbO^. 

The  observation  of  tunable,  stimulated  optical  and  infra¬ 
red  outputs  in  scattering  from  the  248  cm  ^  mode  of  LiNbO^, 
without  the  use  of  an  external  resonator. 

(5)  The  measurement  of  spontaneous  scattering  from  the  248  cm-1 
mode,  whose  frequency  tuning  accurately  confirmed  the  calcu¬ 
lated  dispersion  characteristic,  and  whose  magnitude  varia¬ 
tion  with  frequency  agreed  closely  with  that  predicted  by 
the  developed  theory. 

(4)  The  construction  of  a  spontaneous  scattering  test  set,  con¬ 
sisting  of  an  argon  laser  and  a  phase-sensitive  detection 
system. 

Chapter  II  of  this  report  contains  a  discussion  of  vibrational 
modes  in  crystals  and  the  derivation  of  the  polariton  dispersion  curve 
for  infrared  active  media.  A  brief  discussion  of  phase  matching  on 
the  polariton  modes  is  followed  by  a  presentation  of  a  simple,  physical 
picture  of  the  Raman  process.  Finally,  a  qualitative  comparison  of 
spontaneous  and  stimulated  scattering  is  presented. 
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Chapter  III  contains  the  theoretical  analysis  and  calculations. 

First,  a  classical  analysis  of  the  stimulated  process  is  given.  Using 

2k  25 

an  energy  density  approach,  ’  it  closely  follows  the  method  of 
21 

Henry  and  Garrett,  extending  their  results  to  a  material  with  more 
than  one  vibrational  mode  and  removing  the  "high-loss"  restriction 
made  in  their  paper.  Then,  a  quantum  mechanical  calculation,  using 
rate  equations,  is  performed,  yielding  the  spontaneous  scattering 
cross  sections.  These  results  are  then  related  to  the  stimulated 
gain  coefficient.  This  gives  the  expressions  needed  for  the  calcula¬ 
tions  of  the  nonlinear  constants  of  a  material. 

In  Chapter  IV,  the  results  of  the  previous  chapter  are  applied  to 
LiNbO,.  After  the  calculation  of  the  nonlinear  constants  for  the 

A  -symmetry  modes,  the  stimulated  gain  and  idler  power  are  determined 

^  -1  -1 
for  the  two  major  modes  (248  cm  and  628  cm  ).  The  spontaneous 

scattering  parameters  (for  both  signal  and  idler)  for  these  two  modes 

are  also  calculated. 

In  Chapter  V,  three  sets  of  experiments  are  reported.  First, 
the  argon  laser  and  phase- sensitive  detection  system  are  described, 
along  with  the  study  of  the  idler  power  generated  in  spontaneous  scat¬ 
tering  from  the  628  cm  ^  mode  in  LiNbOy  Then,  the  observation  of 
tunable,  stimulated  optical  and  infrared  emission  in  scattering  from 
the  248  cm  1  mode  is  discussed  and  compared  to  theory.  Finally,  spon¬ 
taneous  scattering  measurements,  which  agreed  with  the  results  of  the 
developed  theory,  are  described. 

Chapter  VI  summarizes  the  highlights  of  this  work.  The  calcula¬ 
tions  performed  in  this  work  are  done  in  esu  units.  The  choice  of  this 
system  was  made  because  the  two  major  references  for  this  work  use 
esu  units.  The  comparison  between  the  theory  developed  in  this  work 
and  that  of  these  two  references  would  have  been  more  difficult  if  a 
different,  system  had  been  used.  For  those  interested  in  converting 

the  results  to  MKS  units,  the  author  re  'oramends  the  use  of  Table  4 

27 

in  the  Apper.d:x  ->f  Jackson. 
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Extensive  use  is  made  of  the  designation  of  frequency,  energy 
and  propagation  constant  in  terms  of  wave  numbers  (cm  ^) .  Conversion 
of  these  quantities  in  cm  ^  to  other  units  may  be  accomplished  with  the 
aid  of  the  Nomograph  of  Electromagnetic  Conversions  of  Pantell  and 
Puthoff1+2(p.  36l). 

The  labeling  of  the  frequencies  of  the  three  waves  is  done  with 
the  subscripts  p  ,  s  ,  and  i  for  the  pump,  signal,  and  idler, 
respectively. 
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CHAPTER  II 


TECHNICAL  BACKGROUND 


A.  vr;  A ! ONAt ,  MODES  IN  CRYSTALS 


Shoo  ‘his  work  is  concerned  with  light  scattering  from  crystals, 
wo  begin  with,  a  brief  discussion  of  the  vibrational  modes  in  crystals. 

train  rerul  *.  of  this  section  will  be  the  calculation  of  the  dis¬ 
tension  curve  for  crystals  with  infrared-active  modes. 

The  solution  h.r  the  propagation  of  a  transverse  vibrational 
wave  through  a  diatomic  lattice  yields  a  dispersion  characteristic 

pQ 

with  two  branches/'  This  is  shown  schematically  in  Fig.  2.1.  The 
lower  branch,  which  intercepts  the  origin,  is  called  the  acoustic 
branch,  since  the  frequencies  associated  with  it  are  approximately 
thi  se  of  sound  waves .  The  upper  branch  is  called  the  optical  branch 
for  the  analogous  reason.  If  the  results  for  the  diatomic  lattice 
are  generalized  to  a  material  with  p  atoms  in  a  primitive  cell,  the 
dispersion  curve  for  transverse  vibrations  in  one  direction  is  composed 

of  one  acoustic  branch  and  (p-l)  optical  branchnes. 

20 

It  can  also  be  shown  that,  on  the  optical  branch,  alternate 
planes  of  adorns  vibrate  against  each  other,  while  the  center  of  mass 
of  th*‘  primitive  cell  remains  fixed.  If  the  atoms  in  alternate  planes 
carry  opposite  charges,  then  vibrations  of  this  type  are  capable  of 
at -sort  ;ng  ar.d  radiating  electromagnetic  energy  at  the  vibration 
frequency . 

The  presence  of  infrared  activity  alters  the  dispersion  curve  of 
tie  rr,a*eria!  in  the  region  of  the  vibration  frequency  due  to  coupling 
*.<  t he  et«vctr  magnetic  wave.  This  will  now  be  derived  since  this 
ohio.  r.  introduces  some  of  the  concepts  and  quantities  which  will 
be  us-  *  >  y‘er.s  .vejy  later. 
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FIG.  2 


K 


7T 

2  a 


1— Dispersion  curve  for  transverse  vibrations  of  a  diatomic 
lattice.  The  lover  curve  is  the  acoustic  branch;  the 
upper  one  is  the  optical  branch.  A  sketch  of  the  motion 
of  the  atoms  during  vibration  is  shown  next  to  each  branch. 
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We  start  with  the  coupled-wave  equation: 

1  S^eM  bir  o^P(co) 

7  x  [7  x  E(ui)]  f  —  eM  - 2  - 2  - -  ’  (2<:L) 

c  St  c  St 

where  S(-u)  and  ?M  are  the  electric  field  and  polarization  oscil¬ 
lating  at.  frequency  id  ,  c  is  the  velocity  of  light  and  is  the 

dielectric  constant  on  the  high  frequency  side  of  the  vibration.  For 
simplicity,  the  problem  for  a  material  with  a  single  mode  is  solved 
and  the  results  are  then  extended  to  the  multimode  case. 

The  polarization  results  from  mechanical  vibration,  which,  in 
turn,  is  related  to  the  electric  field  through  the  equation  for  the 
driven  harmonic  oscillator;  thus, 

PM  =  NeQM  ,  (2.2) 

where  e  is  the  charge  assoicated  with  the  transverse  vibration; 

N  is  the  number  of  primitive  cells  per  unit  volume;  Q  is  the  non¬ 
equilibrium  separation  of  the  vibrating  particles  within  the  primitive 
cell ,  and 


p[Q(u>)  +  TQM  +  uuQM  ]  =  eE(co)  , 


(2.3) 


where  p  is  the  reduced  mass  of  the  vibrating  particles;  is 

the  ;  esor.ar.t  frequency  of  the  vibration;  and  pf  is  the  damping  eo- 
effici-r.t . 

^  Now,  with  the  assumption  that  E(u>)  ,  P(u>)  and  Q(u>)  vary  as 

^i(k  r-tl)  ^  E qs .  (2.i)  and  (2.3)  become  respectively, 


?  *  \ 

-k  +  ~2  D 

c  / 


pM 


(2.M 


ar.c 


where 


mD(u>)Q(^)  »  eE(uj) 


DM  s 


?  2 

i'o  -  -  ?r<» 


(2.5) 


(2.b) 
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Solving  for  Q(cu)  and  substituting  into  Eq.  (2.2),  we  obtain 

2 


P(«i)  = 


Ne 


HD(£d) 


E(a>) 


(2.7) 


Substitution  of  Eq.  (2.7)  into  Eq.  ( 2.k )  gives  the  result: 


2  to  /  U«Ne£ 

k  =  -3  k  +  — 

c  \  pD(cn) 


And  finally  the  definition  of 


puts  the  result  in  the  form: 


to  u 

~ o  ,  €  + - 

c  *  D(o>) 


U) 


e(u>)  | 


(2.8) 


This  is  sketched  in  Fig.  2.2  with  the  real  parr  and  imaginary 
part  of  k  plotted  versus  co  ,  where  k  *  k'  +  ik"  .  The  real 
part  of  k  ,  k'  ,  determines  the  propagation  constant,  while  the 
imaginary  part,  k"  ,  specifies  the  absorption. 

A  second,  perhapamore  physical  way  of  visualizing  the  coupling 
effect  is  shown  in  Fig.  2.3  for  the  lossless  case.  Dashed  lines 
represent  the  separate  propagation  curves  for  the  electromagnetic 
and  vibrational  waves.  The  slope  of  the  electromagnetic  wave  dis¬ 
persion  curve  is  constant  and  inversely  proportional  to  the  index  of 
refraction.  For  simplicity,  the  lattice  optical  branch  is  drawn  as 
a  straight  line.  "Far"  from  the  resonant  frequency  u>0  (below 
point  A)  the  dispersion  curve  follows  that  of  the  pure  electromagnetic 
wave  and  the  curve  is  photon-like.  In  this  region,  the  magnitude  of 
the  vibrational  displacement  is  negligible  and  most  of  the  energy  of 
the  coupled  system  is  electromagnetic.  As  the  resonant  frequency  is 
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FIG.  2. 2- -Frequency  plotted  against  the  propagation  constant,  k' 
and  the  electric  field  absorption  constant,  k"  ,  for  an 
infrared  active  material  with  a  single  mode.  The  right 
hand  figure  is  the  dispersion  curve. 


approached,  (between  points  A  and  B),  coupling  between  the  electro¬ 
magnetic  and  vibrational  wave  increases  and  the  curve  bends  sharply. 

In  this  region,  the  energy  of  the  system  is  divided  nearly  equally 
between  mechanical  and  electromagnetic  energy.  We  will  call  this 
the  coupled  region  of  the  curve.  Beyond  point  B  ,  the  curve  is 
phonon- like  and  most  of  the  energy  is  mechanical.  In  the  lossless 
case,  k  ->  00  as  ui 

In  the  frequency  limit  oi  »  to  ,  the  electromechanical  coupling 
is  negligible  and  the  dispersion  curve  is  that  of  the  electromagnetic 
wave.  As  is  reduced  from  this  limit,  the  coupling  increases  and 
the  curve  flattens  out  again,  becoming  phonon-like.  The  solution  for 
k  =  0  is  easily  obtained  by  setting  e(o>)  =0  in  Eq.  (2.8)  and  it 
is  found  that 


“L 


/^DC  * 

< 


ox. 


> 


where 


'DC 


e 

00 


(2.9) 


2Q 

It  can  be  shown  by  considering  longitudinal  wave  solutions  to 
Eq.  (2.1)  that  the  longitudinal  mode  frequency  is  also  determined  by 
k  =  0  and  therefore 


^0  m‘«L 


(2.10) 


This  is  the  familiar  Lydanne-Sachs-Teller  relationship  for  the 

lossless  case.  The  equivalent  relationship  for  a  lossy  material  has 
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been  derived  by  Barker. 

The  extension  to  a  material  with  more  than  one  vibration  frequency 
is  straightforward.  In  this  case, 


P(cn) 


N.e2 

JLJL 


i 

L 


.1  MjDjWj 


E(o>) 
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and 


(2.11) 


iT  .cd 

Thus,  the  dispersion  curve  for  such  a  material  has  a  "bump"  at 
each  vibration  frequency  ,  the  magnitude  of  the  bump  being 

proportional  to  for  that  mode. 

The  division  o?  the  dispersion  curve  into  three  regions  as 
shown  in  Fig.  2.3,  will  take  on  additional  significance  in  Chapter  III, 
where  the  scattering  theory  is  developed.  There,  the  photon-like, 
coupled,  and  phonon-like  portions  are  relabeled  as  the  parametric, 
vibrational  and  Raman  regions.  This  is  done  to  point  out  the  regions 
in  which  different  terms  in  the  nonlinear  polarizations  play  a 
predominant  role. 


where 


CD 


ft. 


e  + 


l  -bL 


J  D  .  (<u)  * 
3 


— 2  e(“) 


4itN 

_ 


po 


>  S  (u>) 
J 


2  2 
=  O>0  -  CD  - 


B.  THE  SCATTERING  PROBLEM-GENERAL  CONSIDERATIONS 


The  two  conditions  which  must  be  satisfied  in  a  scattering  pro¬ 
cess  are  conservation  of  energy  and  momentum.  We  consider  specifically 
the  three-wave  process  described  in  Chapter  I,  where  a  pump  wave  is 
scattered  from  an  optical  mode  in  a  crystal,  generating  signal  and 


it  nota- 
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Barker  and  Barker  and  Loudon  use  a  slightly  u. 

2  2 

tion.  They  define  S  ,cd.  s  JJ  ,  where  S.  is  called  the  oscilla- 

0  °j  Pj  J 

tor  strength  for  the  mode.  This  form  yields  the  result 


eDC  “  €°° 


+  5 


(2.12) 
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idler  waves.  The  appropriate  equations  for  the  interaction  are 


where 


OJ  =02 

P  s 

+  U2. 

1 

(2.13a) 

t'  =  t' 

P  s 

+  , 

(2.13b) 

k'  =  23. 

c 

. 

The  vector  diagram  for  Eq.  (2.13b)  is  shown  in  Fig.  2.4,  where 

0  is  the  angle  between  It  and  It  and  varies  from  0°  for  the  col-> 

P  0  s 

linear  forward  scattering  to  180  for  collinear  backward  scattering. 

This  work  is  main^„  concerned  with  scattering  for  0  <  10°  which  is 

generally  referred  to  as  near-forward  scattering. 

Solution  of  Eq.  (2.13b)  together  with  Eq.  (2.13a)  yields  a 

family  of  k^  vs  V .  curves  for  different  values  of  0  .  These 

curves  are  obtained  from  Eq.  (2.14),  with  the  substitution 

V .  =  v  -  v  •  The  intersections  of  these  curves  with  the  dis- 
i  p  s 

persion  curve  indicate  the  frequency  and  wave  vector  of  the  idler 
which  will  take  part  in  the  scattering  process  for  each  angle. 

This  is  shown  schematically  in  Fig.  2.5. 


C.  THE  RAMAN  EFFECT  -  A  SIMPLE  PHYSICAL  PICTURE 
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In  this  section,  a  simple  physical  picture  of  the  interaction 
of  light  with  a  vibrational  mode  is  given. 

The  pump  is  at  frequency  cjo  which  is  assumed  to  be  far  above 
the  highest  vibrational  resonance  of  the  molecular  ions  (i.e., 

0)  »  u>A  )  .  To  zero  order,  therefore,  the  ions  do  not  respond  to 

p  Uj  • 

the  applied  field  and  only  the  electrons  contribute  to  the  polariza¬ 
bility.  Neglecting  the  dispersive  effect  of  the  "tails"  of  the  elec¬ 
tronic  resonances,  the  polarizability,  a  ,  will  be  independent  of 
frequency.  The  polarization,  given  by 


P  *  «E  = 


otE  cos  as  t 
P  P 


(2.15) 
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i\|  *  r\p“ '*5 

K’j2  *Kp  +  K*  -  2Kp  Kg  cos  9 


*47T2(j/pT}p-i/8  i;8)2+8T^z/pi/8  ?78(I-cor0)  (2.14) 

WHERE  V  IS  THE  FREQUENCY  IN  cm  1 
AND  T)  IS  THE  REFRACTIVE  INDEX 


FIG.  2. 4 — Momentum  conservation  triangle.  The  angle  9  is  the  phase¬ 
matching  angle.  The  solutions  of  (2.13b),  together  with 
the  energy  conservation  condition  (v^  =  vp  -  vg)  ,  are  a 
family  of  curves  of  vs  Vj^  .  These  curves  are  referred 
to  as  phase-matching  lines. 
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FIG.  2. ^--Schematic  representation  of  the  phase-matching  lines, 
superimposed  on  the  dispersion  curve  for  a  two-mode 
material.  The  points  of  intersection  indicate  the 
idler  frequencies  participating  in  the  scattering 
process.  For  G  >  ,  the  phase-matching  lines  do 

not  intersect  the  dispersion  curve,  and  the  idler 
is  of  phonon  nature  in  this  region.  The  idler  frequency 
asymptotically  approaches  the  resonant  vibration 
frequency  for  large  •)  .  This  has  been  verified 

experimentally. » 32 
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will  therefore  contain  only  a  term  at  a> 

P 

There  is  some  energy  transfer,  however,  between  the  electron 
cloud  and  the  ions,  causing  small  ionic  vibrations.  In  this  case, 
the  electronic  polarizability  is  a  function  of  the  coord' nates  of  the 
ions.  If  the  polarizability  is  written  as  an  expansion  in  the  normal 
modes  of  vibration  about  the  equilibrium  position,  we  get 


da 

a  =  aA  +  Y.  -  Q.,  +  higher  order  terms  ,  (2.l6) 


where 


Q.  =  COS  Oj.  t 

J  0  lj 

Now,  substitution  of  Eq.  (2.l6)  into  Eq.  (2.15)  yields  (to 
first  order) 

_  da 

P  *  a.E  +  L  -  Q.E 

0  ji«,  1 


^  da  _  _ 

=  a_E  cos  a>  t  +  i,  -  Q.E  cos  (a>  -a>.  )t 

Op  P  H  »  *3 

da 

+  T  -  Q.E  cos  (a)  )t 

3  1  p  p 


(2.17) 


Thus  polarizations  are  generated  at  the  sum  and  difference 

frequencies.  Historically,  the  line  at  the  difference  frequency  (u>  ) 

is  called  the  "Stokes"  output  and  the  line  at  the  sum  frequency  (u/) 

8 

is  called  the  "anti-Stokes"  output.  In  spontaneous  scattering  and 
under  most  conditions  in  stimulated  scattering,  the  Stokes  output  is 
stronger  than  the  corresponding  anti-Stokes  output. 

In  this  discussion,  u>^  has  been  taken  to  be  the  resonant  frequency 
of  a  normal  mode  of  vibration.  As  described  previously,  in  infrared- 
active  media,  the  vibrations  can  be  driven  off  resonance.  Hence,  <1^ 
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can  be  generalized  to  represent  any  frequency  on  the  dispersion 
curve  for  the  material.  As  long  as  the  wave  vector  and  frequency 
conservation  conditions  are  satisfied,  tunable  outputs  at 
and  -jo.  can  be  generated. 


D.  QUALITATIVE  COMPARISON  OF  STIMULATED  AND  SPONTANEOUS  SCATTERING 


While  both  stimulated  and  spontaneous  scattering  arise  from  the 
same  physical  process  (described  briefly  in  the  previous  section), 
there  are  differences  between  them.  Qualitatively,  some  of  the 
important  differences  are  as  follows: 

* 

a.  Stimulated  scattering  is  a  threshold  process,  requiring 

8  ,  2 

power  densities  on  the  order  of  10  W/cm  to  initiate  the 
process.  Power  densities  of  this  magnitude  are  available 
only  from  Q-switched  lasers.  It  is  for  this  reason  that 
stimulated  scattering  was  r.ot  observed  until  after  the 
development  of  the  Q-switched  laser. 

b.  Signals  generated  in  stimulated  scattering  aure  of  narrow 
spectral  width  and  are  well  collimated.  In  contrast,  spon¬ 
taneous  scattering  outputs  are  scattered  into  all  directions 
and  possess  the  spectral  width  of  the  scattering  vibrational 
mode . 

c.  The  conversion  efficiency  (generated  signal  power  divided  by 

input  pump  power)  of  spontaneous  scattering  is  about  10  , 

while  the  stimulated  scattering  conversion  efficiency  is 
typically  1-5$  and  theoretically  can  approach  10$. 

These  and  other  properties  of  the  two  processes  are  presented  for 
comparison  in  Fig.  2.6.  Many  of  these  will  be  discussed  quantitatively 
in  later  chapters. 


<1 

For  the  case  we  are  considering,  the  terms  stimulated  and 
threshold  are  defined  as  follows:  the  signal  emission  is  stimulated  when 
the  exponential  growth  term  predominates  over  the  linear  growth  term. 

The  power  density  which  produces  this  is  the  threshold  power  density. 
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FIG.  2.6--A  qualitative  comparison  of  some  of  the  properties  oi 
stimulated  and  spontaneous  scattering. 
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CHAPTER  III 


THEORETICAL  ANALYSIS 


A.  INTRODUCTION 

This  chapter  is  devoted  to  a  theoretical  analysis  of  stimulated 
and  spontaneous  scattering.  First,  a  classical  analysis  of  the  stimu¬ 
lated  scattering  problem  is  carried  out.  Using  the  energy  density 

2i>  26  21 

approach,  *  the  analysis  follows  the  method  of  Henry  and  Garrett. 

The  significant  contribution  of  the  present  work  is  to  extend  their 
results  to  a  multimode  material  and  to  derive  a  gain  expression  which 
is  not  restricted,  as  theirs  is,  to  either  the  collinear  case  or  the 
"high-loss"  case. 

A  quantum  mechanical  calculation,  using  rate  equations,  follows. 

When  applied  to  the  stimulated  process,  the  results  of  this  analysis 
agree  with  the  classically  derived  expressions  as  expected.  In 
addition,  the  scattering  parameters  of  the  spontaneous  process,  which 
can  not  be  arrived  at  classically,  are  derived.  Finally,  the  procedures 
for  computing  the  material  nonlinear  coefficients  are  determined. 

B.  STIMULATED  SCATTERING  -  CLASSICAL  ANALYSIS 

A  classical  analysis  of  stimulated  scattering  is  presented  in 

I'. 

this  section.  This  problem  has  been  discussed  previously  by  Loudon,  ' 

18  19  21 

Butcher,  Loudon,  and  McLean,  shen,  Henry  and  Garrett  and  in  an 

qU 

unpublished  work  by  Puthoff.  The  calculations  presented  in  this  sectiqn 
are  based  upon  the  work  of  the  latter  two  references.  In  particular,  the 
approach  of  Henry  and  Garrett,  hereafter  called  HG  ,  will  be  followed. 

We  begin  by  assuming  plane  waves  for  the  pump,  signal  and  idler 
electromagnetic  fields  (E)  and  for  the  lattice  displacement  field  (Q). 
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Thus 


i(k  -r  -  as bt)  +  y  -r 

Ee  '  e;<V  *  c-°-  '  As  e  +  c'°-  ' 

(3-1) 

where  (i  =  p  ,  s,  i  for  the  pump,  signal  and  idler  fields  respectively, 
y  is  the  growth  (or  decay)  constant,  and  c.c.  denote.-  complex  conjugate; 
and 

iO^  r  -  o^t)  +  7i* r 

=  Q(o)i)  +  c.c.  =  Qj  +  c.c. 

(3*2) 

Before  proceeding,  several  comments  concerning  the  definition  of 
the  plane-wave  fields  are  ir.  order.  A  possible  source  of  confusion  is 

the  definition  of  the  complex  amplitudes  without  the  usual  factor  of  j. 

35 

This  convention,  first  introduced  by  Pershan,  is  used  here  because 
HI  use  it.  This  leads  to 

Re  (E^)  =  2  Ap  cos  (kg  *  r  -  u^t)  , 


or  a  real  amplitude  of  2A^  .  This  has  two  consequences  of  which  one 

has  to  he  aware  in  order  to  avoid  the  annoying  "factor  of  two"  diffi¬ 
culties.  First,  the  intensity  of  the  plane  wave  defined  in  is 


h 


1^1 


(30) 


(The  more  common  expression  uas  a  factor  of  8  in  the  denominator. ) 

Secondly,  as  discusaed  by  Bloembergen,  1  this  definition  leads  to 
a  factor  of  increase  in  the  lowest  order  nonlinear  susceptibilities 
and  additional  factors  of  2  increase  in  higher  order  ones.  Thus,  in  the 
body  of  this  work,  when  results  are  compared  to  those  in  the  literature, 
(and  to  Putboff,^4  who  uses  the  more  widespread  convention),  the 
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appropriate  factors  must  be  introduced.  This  will  be  noted  each  time 
it  is  done  to  minimize  confusion  for  the  reader. 

To  continue  with  the  analysis,  we  note  that  the  three  basic 
equations  used  to  describe  the  problem  have  already  been  referred  to 
in  Chapter  II.  They  are  the  coupled  wave  equation  (2.1 j,  (written 
for  the  three  fields  of  interest),  the  polarization  equation  (2.2), 
and  the  driven  oscillator  equation  (2.3).  The  significant  difference 
between  the  analysis  of  that  chapter  and  the  analysis  here  is  that 
in  Chapter  II  attention  was  confined  to  the  solution  of  the  linear  problem, 
by  using  only  the  linear  polarization  term  in  the  coupled-wave  equation 
and  the  linear  f~_ ce  term  in  the  oscillator  equation.  In  this 
chapter,  the  nonlinear  driving  terms  will  be  included. 

The  coupled  wave  equation  is  written  in  (2.1).  In  general,  with 
fields  polarized  in  more  than  one  direction,  the  tensor  character  of 
the  dielectric  constant  must  be  included.  In  this  work,  however,  we 
are  concerned  with  the  situation  where  all  fields  are  polarized  in  the 
same  direction,  (along  the  c-axis  of  the  crystal).  Therefore,  the 
tensor  notation  is  not  necessary.  In  addition,  it  should  be  pointed 
out  that  the  polarization  in  (2.1)  includes  both  linear  and  nonlinear 
terms  where  applicable. 

A  generalization  of  the  driven  oscillator  equation  (2,3)  yields 
Dd  (».)  Q.(a,.)  .  Fjfcj)  »  e3  !(»,)  +  F®1^)  , 

(3-4.) 

f  y. 

where  the  equation  has  been  written  for  the  j  vibrational  mode  of 
a  multimode  material  and  the  total  force  has  been  divided  into  two 
parts;  the  linear  part,  (which  is  the  same  as  before),  and  a  nonlinear 
term  arising  from  the  scattering  process. 

There  are  several  approaches  to  the  calculation  of  the  nonlinear 
driving  terms for  (2.1)  and  (3.U).  One  approach,  developed  by  Puthoff,^ 
is  to  derive  the  equations  of  motion  for  the  dipole  moments  and 
polarizability  from  quantum  mechanical  considerations  using  the  density 
matrix  approach.  These  equations  are  then  solved  to  yield  the 
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polarizations  in  terms  of  susceptibilities  which  are,  in  turn,  defined 
in  terms  of  the  matrix  elements  of  the  electric  dipole  moments  and 
polar i zab ilit ies . 

A  second  method,  which  is  used  by  HG  and  which  will  be  used  here, 
assumes  the  existence  of  a  phenomenological  energy  density  function,  U  , 

from  which  the  nonlinear  driving  terms  can  be  derived  directly.  This 

2  5 

approach  was  first  suggested  by  Kleinman.  A  detailed  discussion  of 
the  energy  density  approach  is  presented  in  Appendix  A. 

In  the  problem  under  consideration,  there  are  electromagnetic 
fields  present  at  m  ,  05g  and  ,  and  a  lattice  displacement  at  , 

with  contributions  to  the  lattice  displacement  from  j  modes.  The 
appropriate  form  of  the  energy  density  fuction  is  shown  in  Appendix  A 
to  be 


UtQ^),  E(oy.),  E(ods)  ,  E(o>p)3 

=  -  [dg  E(<up)  E(o>s)*  E(o>i)*  +  E(o)p)  E(u5g)*  1  N j  ^  +  e-C. 

(3-5) 


The  quantities  d^,  and  d^  are  the  nonlinear  coefficients  for  the 
material.  It  has  been  shown^  that  the  taking  of  dg  and  d^  to  be 
real  and  constant  is  consistent  with  experimental  results.  Also,  it 
shall  be  demonstrated  later  that  d^  is  related  to  the  more  familiar 
nonlinear  susceptibility.  The  quantities  d^,  and  d^  are  third  rank 
tensors,  but,  «ts  explained  before,  the  tensor  notation  is  not  needed 
to  describe  the  process  of  interest  for  this  work. 

As  indicated  in  Appendix  A,  the  nonlinear  driving  terms  are 
derived  from  the  energy  density  function  as  follows: 


(3.6a) 
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and 


rfW  - 


J  1 


su 


n3 


(3.6b)’ 


Taking  the  appropriate  derivatives  and  substituting  into  (2.1) 
and  (3.4),  we  obtain 


[,2  “i 

!v  +  ~  e« 
c  i. 


\ 


4rt<u. 


E(o>.)  = 


cL,  E(o)p)  E (a>s )*  +  Y.  Nj  sj  K) 

J 


(3.7a) 


(f 


a) 


+  1e  E(o)s) 
\  c  s  / 


iV 

\  c2  V  ~'~P 


4jtm2  P 

S 


E(cn  )  =  - 


4«cd' 

"2 


[ds  ES}  E(a,i}*  +  ES)^  Nj  ^ 


(3.7b) 


2  _ 


j^E  E(u>s}  EK)  +  EK}  T,  Nj  dQj  W 


(3.7c) 


Uj  (ax) (ax )  =  e^  E(ax)  +  E(o)p)E(a»s)*  .  (3«7d) 


The  solution  of  (3.7  a-d),  subject  to  the  appropriate  boundary 
conditions,  yields  a  complete  description  of  the  process. 


Note  that  (3)  of  HG  is  incorrect.  A  simple  dimensional  analysis 
will  show  that  the  equation  above  is  the  correct  one.  It  is  appropriate 
to  point  out  at  this  time  that  the  middle  equation  of  (2)  of  HG  is  also 
a  misprint  since  the  right-hand  side  should  contain  the  total  polarization 
(nonlinear  and  linear  terms).  Both  of  these  misprints  have  been  eliminated 
in  (7)  of  their  paper  and  hence  their  results  are  correct. 


-  26  - 


In  practice,  however,  it  is  impossible  to  obtain  a  complete 
solution  in  closed  form  without  making  some  assumptions.  The  first 
assumption  is  tc  neglect  pump  depletion.  This  is  equivalent  to  the 
statement  that  the  polarization  terms  on  the  RHS  of  (3* 7c)  are  small 
and  do  not  significantly  alter  the  intensity  of  the  pump  field  from 
its  initial  value  at  the  input  of  the  crystal.  Therefore,  (3- 7c) 


is  drqpped,  A  is  treated  as  a  constant  and 


y  =0 
P 


Depletion 


of  the  pump  does  occur  under  certain  conditions  in  the  experiments 
to  be  described  later  and  therefore  this  assumption  is  not  completely 
valid  in  that  case. 

To  gain  more  physical  insight  into  the  problem,  we  solve  for  Q.  (a>. ) 

J 

in  (3. 7d)  and  substitute  this  result  into  (3.7  a,b).  This  yields 


(<V*  - 


!>l  Etu,)* 

2j _ _ 


ap,  dS,  E(o) I  )*  K(a>  ) 

-J - i - E - —  (3.8a) 


2 

2  * 

^+~€i 

c 


2 

“i 


®K)  -  -T 


-  p  / 

ap  ^  1 

E(cd  )*  E(tDs)  +  E((i)J*  e((Ds)  . ( 


j  D  (  a,i)*J 


(3.8b) 


2  1 


co 
,  s 
+  ~o  e 

C-  co 
C  S 


E(cos) 


2  ' 

S.V 


d£  E(co  )  E(co.f  ^  ~JL±  jE(co  )  ECco.f 


j  (oi) 


*  J 


O  'O 

ap,  % 

+  |E(co  )|2  eV)  )  i 

p  s  L-j  ,  .*  , 
i  Dj(  wi)  J 


(3»8c) 
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To  get  these  equations,  we  have: 

a)  introduced  the  simplifying  notation  of  HG: 


,  dg  _  4itdg 


Si  a  4”  HJ  e3  **3 


b)  written 


..♦I 


c)  taken  the  complex  conjugates  of  (3*?a)  and  (3-7d)  . 

It  is  instructive  to  consider  the  three  driving  terms  on  the  RHS 
of  (3.8c).  The  first  term,  d^  E(a^)  E^)  ,  is  a  typical  polarization 
associated  with  parametric  processes.  It  depends  only  on  the  product 
of  two  electromagnetic  fields  and  exists  independently  of  any  vibrations. 
It  will  be  largest  in  the  region  where  E(cd^)  is  largest.  As  was 
pointed  out  earlier,  this  will  be  in  the  photon-like  portion  of  the 
dispersion  curve  where  most  of  the  energy  of  the  mixed- system  is 
electromagnetic.  For  this  reason,  the  photon-like  portion  of  the 
dispersion  curve,  (between  the  origin  and  point  A  of  Fig.  2.3),  is 
also  called  the  parametric  region. 


The  last  term. 


V’  P1  c 

L  ~L~i 

5  DjUp 


is  the  normal  Raman  polarization  term.  It  originates  from  the 
product  of  the  pump  field  and  the  lattice  vibration  displacement 
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and  in  no  way  depends  upon  the  existence  of  E(an)  ,  (i.e.,  it  can  exist 
in  materials  which  are  not  infrared  active).  This  term  will  predominate 
where  Q(ox)  is  large,  which  is  in  the  phonon-like  portion  of  the  dis¬ 
persion  curve.  In  addition,  the  resonant  denominator  maximizes  this 
term  at  ;  therefore,  the  "pure"  Raman  gain  is  largest  at 

cn.  =  o>p  -  ^  and  this  is  the  frequency  at  which  "pure"  Raman  scat¬ 
tering  is  always  observed.  Hence  the  phonon-like  portion  of  the  curve, 
(between  point  B  and  the  vibration  frequency),  is  called  the  Raman 


region. 


The  middle  term. 


4  dQ 

Q. 


J  DjH)* 


E(i»p)  Efo,;' 


depends  upon  the  existence  of  both  a  lattice  vibration  and  an  electric 
field  at  ox  .  The  influence  of  this  term  is  felt  most  strongly  in 
the  middle  (or  coupled)  region  of  the  dispersion  curve  due  to  the 
counterbalancing  effects  of  E(cui)  and  D^(cd^)  .  This  term  has  been 

called  the  vibrational  term  and  the  corresponding  region  of  the  dis¬ 
persion  curve  (between  points  A  and  B)  has  been  labeled  the  vibra¬ 
tional  region 

At  this  point  we  compare  the  coupled-wave  equations  for  cu.  and 

oO  1 

o)  [(3.8b)  and  (3.8c)]  to  those  given  by  Puthoff.  A  term-by-term 
s 

comparison,  taking  into  account  the  comments  after  (3-2),  reveals  the 
following  equivalences : 


Present  Work 


Puthoff 


d'2  a2 


(3.9a) 


(3- 9b) 


DjK) 


(3-9c) 
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where  Xy  ,  X^  and  Xp  have  been  identified  by  Puthoff  as  the 
vibrational,  Raman  and  parametric  susceptibilities  respectively. 
(Note  that  these  definitions  do  not  correspond  to  our  usage  of  the 

ol| 

same  terms).  Referring  to  (4.5)  and  (4.6)  of  Puthoff,  we  find 
that  the  vibrational  and  Raman  susceptibilities  can  be  written  as: 


rl  “0,  V,  <«b  > 

Xv_(o).)  =  - * - * - 3- 


rj  “o,  %R.  H>  ^ 

v  (a,  )  =  - * - L - , 

^  K) 

where  x."  (con  )  is  the  value  of  the  imaginary  part  of  the  susceptibility 

j  u . 

evaluated  on  resonance  (ax  =  o)q  )  .  Using  these  expressions  and 

(3*9a>b),  we  establish  the  following  relationships: 


■  *  ra  <»0j> 


\  "pj  •  16,,ra  "0^ 


(3.10a) 

(3.10b) 


The  latter,  in  particular,  will  be  useful  in  later  calculations. 

Returning  to  the  solution,  the  plane-wave  expressions  for 
E^)  ,  E(cl>s)  and  Q(ax)  [(3.1,  3.2)]  are  substituted  into  (3.7). 
Exponential  balance  of  both  sides  of  the  equations  requires: 

k?  =  k.  +  (3- Ua) 

7i  =  7S  =  7  .  (3* Hb) 

The  simple  result  of  (3.11b)  is  due  to  the  assumption  of  no 
pump  depletion  (7  =  0).  Cancellation  of  the  equal  exponentials 
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CM  OT 


reduces  (3.7)  to: 


For  a  solution  to  exist, 


We  now  define  the  gain  as  the  intensity  change  along  the  direction 
of  propagation,  which  is  very  nearly  the  direction  of  the  Poynting 
vector.  Thus, 

gi  =  2^.7  ,  gg  =  2kg-7  ,  (3.13) 
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where  k^  and  kg  are  unit  vectors  in  the  k^  and  kg  directions, 
respectively. 

Expanding  (k.  +  i?)  and  (k  -  i?)  to  first  order  in  y  , 

X  s 

we  obtain 

(k.  +i yf  =  k^  +  2i  k.-7  =  k^  +  iki  6i 
(kg  -  i7)  -  kg  -  i  kg  gg  . 

Substitution  of  this  into  (3-12),  yields  the  final  equation: 


d^  ,  d^  ,  Qp  )  ,  the  frequencies  and  wave  vectors  of  the  inter¬ 

acting  signals  and  the  intensity  of  the  pump,  there  are  still  two 


unknowns- g.  and  g  -and  only  one  equation.  Thus  an  additional 
X  s 

equation  relating  g,  and  g  is  necessary.  1C  solve  the  problem 

X  s 

by  treating  tvo  limiting  cases : 

1)  High-loss  region:  In  this  region,  ( 3 . 14 )  is  simplified  by 


setting  a  0  .  This  can  be  justified  mathematically  in  regions 

of  high  loss  at  the  idler,  where  k^  gi  can  be  neglected  compared  to 


2 

mi 

-x 

c. 
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in  (3.IU).  In  this  region,  the  calculation  of  g  is  done  separately 

s 

for  two  regions  of  for  a  multimode  material: 

a)  For  outside  the  linewidth  of  each  of  the  j  modes,  an 
analysis  is  given  in  Appendix  B,  which  shows  that: 


and 

b)  For  within  the  linewidth  of  the  mode  with  resonant 
frequency  cu^  ,  the  two  summations  in  (3.14)  are  dominated  by  a 
single  term;  ?.e.. 


dT 

i 


< 


DaK}* 


Thus,  in  this  region,  the  calculation  is  identical  to  the  single-mode 
treatment  of  HG  and  hence: 


-  33  - 


If  03. 


U3Q  is  substituted  into  (3.15b),  the  result  is: 
°a 


8°(V  ' 


03  n 

8  P 


a 


'V8  < 


c  “c  \  ra 
a 


(3.15c) 


Using  the  equivalences  of  (3.10b),  this  is  seen  t'  agree  with 
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the  standard  expression  for  the  on  resonance  Raman  gain. 

Thus,  (3.15a)  and  (3.15b)  are  the  extensions  of  (18)  of  HI 
to  a  multimode  material  in  the  high-loss  region.  As  we  shall  see 
later,  this  assumption  is  not  valid  throughout  the  entire  range  of 
the  experimental  observations. 

2)  Collinear  Interaction:  For  the  collinear  process  (e  =  0), 

g.  =  g  =  g  .  The  result  in  this  case  (which  is  a  straightforward 
1  s 

generalization  of  the  analysis  in  Appendix  B  of  H3)  is 


(3*16) 


where  =  2k^  =  power  absorption  coefficient  at  ;  and 


U3.  U3  p 

1  S-  |A  ,2 


Uc2  \  \  P 


H* 


I 


(3.17) 


where  T,  and  T)  are  the  refractive  indices  at  03.  and  03  , 

is  is 

respectively.  Again,  phase-matching  has  been  assumed. 

It  is  of  interest  to  consider  limiting  cases  of  the  collinear 
result : 

a)  Low-loss  region 

In  this  region,  ot^  <-<  ap  and 


ge  =  2a  .  (3.18a) 

8  P 
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and 


b)  High-loss  region 
In  this  region, 


(3.18b) 
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This  can  be  shown  to  be  equivalent  to  the  previous  high- 
loss  expressions  (3.15a)  and  (3.15b). 

The  collinear  asumption  is  obviously  not  strictly  valid  except 
at  the  point  of  intersection  of  the  0=0  phase-matching  line  and 
the  dispersion  curve  for  the  material.  In  addition,  as  we  shall  point 
out  in  more  detail  later,  the  collinear  process  is  not  possible  for 
the  lowest  frequency  mode  of  most  materials  (including  LiNbO^), 
since  the  intersection  of  the  0=0  phase  matching  line  and  the 
dispersion  curve  near  the  lowest  mode  occurs  at  tu.  =  0  . 

X 

Thus  far,  we  have  derived  expressions  for  the  stimulated  gain 
in  multimode  media  in  both  the  high-loss  approximation  and  in  the 
collinear  case.  To  the  extent  that  these  particular  constraints 
apply,  the  expressions  describe  the  stimulated  scattering  process. 

As  has  been  pointed  out,  however,  neither  of  these  conditions  is  valid 
for  the  entire  range  of  experiments  to  be  described.  The  analysis 
which  follows  leads  to  a  more  general  gain  expression,  which  is  valid 
over  the  entire  region  of  near-forward  scattering. 

The  approach  in  this  analysis  is  to  note  that  if  the  direction 
of  y  can  be  determined,  then  the  desired  relationship  between 
and  gs  can  be  obtained  for  all  phase-matching  angles.  This  is 
accomplished  by  considering  the  boundary  conditions  at  the  surface 
of  the  nonlinear  material.  The  details  of  the  calculation,  first 
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suggested  by  Bloembergen  and  Perehan,  are  given  in  Appendix  C. 

It  is  also  shown  in  this  appendix  that  the  desired  relationship 
between  g^  and  gg  for  small-angle  scattering  and  near-normal 
incidence  of  the  pump  is  given  as 


g j  =  gg  cos  ® 


(3.19) 
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where  ip  is  the  angle  between  the  idler  and  pump  propagation  directions 
(see  Fig.  2.!s ). 

As  shown  in  Appendix  D,  this  leads  to  the  final  result  for  small- 
angle  scattering: 


a. 
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2  cos  cp 
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16a 


?  Y/2 
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af 

i 
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2  2 
where  a '  =  cos  cp  a 


Examining  the  same  limits  as  before,  we  find: 


(3-20) 


Lew  Loss: 


2a' 
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2a 


cos  cp  (cos  cp) 


1/2 


(3.21a) 


High  Loss:  g 


P 


cos  cp 


ha 


a. 


(3.21b) 


Note  that  for  collinear  scattering  (cos  cp  =  l),  these  results 
agree  with  those  of  the  previous  analysis.  Also  of  interest  is  the 
fact  that  she  high  loss  result  again  agrees  with  the  previous  expressions 
(3.I5)  and  (3.18b).  The  reason  for  this  is  that,  in  the  high-loss 
region,  the  result  for  g  is  independent  of  the  value  used  for  g. 
since  the  term  involving  g.  can  be  neglected  in  (3. 14). 

The  results  of  this  section  are  summarized  graphically  in 
Fig.  3-1  where  g  is  sketched  as  a  function  of  to.  using  the 
derived  expressions.  Similar  curves,  using  the  calculated  numerical 
values  for  LiNbO^,  are  presented  in  Chapter  IV. 
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FIG.  3.1--Sketch  of  results  of  stimulated  gain  calc  illations,  showing 

signal  gain  coefficient  (g  )  plotted  against  idler  frequency. 
Curve  (b),  drawn  solid,  is  a  plot  of  (3.20)  ana  is  valid  over 
the  entire  mode.  Shown  for  comparison  are  Curves  (a)  and  (c), 
which  are  plots  of  (3-15)  and  (3.16),  respectively.  These  are 
not  valid  over  the  entire  mode.  Note  that  the  three  curves 
coincide  in  the  high-loss  region;  the  reason  for  this  is 
explained  in  the  text. 


C.  QUANTUM  MECHANICAL  ANALYSIS  -  STIMULATED  AND  SPONTANEOUS  SCATTERING 


The  previous  section  contains  a  classical  analysis,  which  is 

valid  for  the  stimulated  process  since  the  quantum  numbers  are  large. 

To  include  spontaneous  effects,  however,  a  quantum  mechanical  analysis 

is  required.  This  is  considered  in  this  section. 

i+2 

We  begin  with  the  rate  equations 


d<n  > 
s 


dt 


=  K(<rn>  +  l)  (<n^>)  (<ns>  +  l)  -  K<hi>  <ng>  (<n^>  +  l) 


(3.22a) 


d<n  >  d<n  >  d<5i.>  <n. > 

_ s_  _ _ p_  _  _ x_  +  _ x_ 


(3.22b) 


dt 


dt 


dt  t 


where  <h>  is  the  number  of  photons  in  a  single  mode,  and  the  sub¬ 
scripts  refer  to  the  frequencies.  Note  that  losses  at  the  pump  and 
signal  have  been  neglected  and  less  at  the  idler  has  been  accounted 
for  by  including  the  term  <tu>/t  .  We  have  also  defined 


where  v^  is  the  wave  velocity  at  and  (X  is  the  absorption 

coefficient  at  00. 

1 

These  equations  are  now  applied  to  stimulated  and  spontaneous 
scattering. 

1.  Stimulated  Scattering 
In  steady  state, 


dt 


v«v 


v-V 
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Therefore,  (3.22b)  becomes 


v  •  7  <ii  >  =  v.-V  <n.>  +  CC.  <n.> 
s  s  1111 


(3*23) 


2.y*y'  t  \ 

Now  assuming  that  <h.>  ,  <h  >  grow  as  e  as  before,  (3*23) 

1  s 

becomes 


2v  .7  <ii  >  =  2v.  •  7  <n. >  +  a.  <h. 

s  '  s  1  '  1  11 


.> 


But 


2v  •  7  =  v  g 

s'  s  °s 


2v.  •  7  =  v.  g. 

1  '  1  Bi 

Therefore,  a  relationship  between  <ni>  and  <hg>  is  derived; 


v  /  g 

,  ^  s  •'  &s  ,  ,  . 

<n.  >  =  —  i  -  <n  > 

1  __  v  _  .  ~  /  s 


v.  \g.  +  a 
1  °i 


(3.24) 


To  calculate  the  idler  power,  we  recall  that 
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<h> 

V 


£L 

ftCOC 


(3.25) 


where  I  is  the  intensity  and  V  is  a  volume  whose  linear  dimension 
L(V  =  It)  is  large  compared  to  a  wavelength,  but  small  compared  to  the 
growth  of  the  wave  (i.e.  L  «  ^)  . 

•  s 

Using  (3.24)  and  (3*25)  yields 


_i _ 1 

P  L 
s  s 


®i 


0) 


[  «■ 

V  g.  +  CC. 
V<51  1 


where  P.  and  P  are  the  generated  idler  and  signal  power  respectively, 

IS  ' 

and  where  31  and  are  "small"  distances  along  the  direction  of 
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propagation  of  cd. 


and  03s  respectively.  But 


L. 

1 


cos  tp 


from  geometric  considerations;  and  therefore,  after  using  (3*19) 
we  get  the  result: 


P. 

1 


“i 


,o> 


g3  cos  cp 


v6s  COS  cp  +  ct± 


(3-?6) 


Examining  the  limiting  case  of  (3.26),  it  is  found  that  in  the  low-loss 
region  (a^  -  0) 


CD. 

P.  =  —  P 
1  s 

ms 


the  Manley-Rowe  limit.  In  the  high-loss  limit, 


P. 

1 


cd.  g  cos  tp 
i  s 


CD_ 


a. 

1 


We  now  calculate  the  signal  gain  in  the  high-loss  limit.  In 
this  limit,  the  proper  assumptions  for  the  photon  numbers  are 


<n  >  ,  <n  >  »  1 
p  ’  s 

»  <n,  > 
1 


Insertion  of  these  assumptions  into  (3.22a)  immediately  yields 


K<n  > 

_ L. 

c/n. 


(3.27) 


Note  that  in  the  high  loss  region  g_  is  proportional  to  pump 

s 

intensity  (a  <n^>),  i'1  agreement  with  the  classically  derived  result. 
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2.  Spontaneous  Scattering 


The  assumptions  for  the  photon  numbers  <h>  for  the  three  waves 
in  spontaneous  scattering  require  a  more  detailed  explanation.  As 
in  the  stimulated  case,  <h^>  >>  >  having  a  value  of  ~  10,+  for  a 

1  watt  pump.  Since  spontaneous  scattering  efficiencies  are  typically 
~  Id'6,  it  may  appear  plausible  to  set  <h  >  ,  <h.>  «  1  .  However, 

<h>  is  composed  of  two  components;  one  due  to  scattering  and  the 
other  arising  from  the  Bose  distribution  (black  t'dy)  which  is  given 

by 


n  = 


.WM  _  ! 


(3-28) 


where  k  is  Boltzmann's  constant  and  T  is  the  temperature  in  K. 
For  room  temperature,  n 


0  ,  since  o)g  is  an  optical 


frequency  and  ho>  »  kT  .  On  the  other  hand,  n.  may  be  equal 

S  1 

to  or  greater  than  1  depending  upon  the  particular  value  of  . 

Therefore,  the  proper  assumptions  for  spontaneous  scattering 


are 


V  >>  1 

"V  =  ^s  +  ns(zs)><:<1 

<ii.>  =  <n.  +n.(z.)>  =  n. 

i  i  i  i  l 


Assuming  the  steady  state  and  inserting  the  assumptions  above, 
we  find  that  (3.22a,b)  become 


cKn  > 
s 

dz_ 


K  (<h.  >  +  1)  <h  >  a 
i  P 


-  V 


> 

_ p_ 

8z 


8<h.> 

v.  - —  +  a.  v.  <n.> 

1  ~dz.  111 


(3.29) 
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To  arrive  at  an  expression  for  the  scattering  cross-section, 
we  solve  (3-29)  for  n^  assuming  <n^>  =  =  constant.  This  gives 


K(n.  +  1  )Ti  z 
_ IP  1 


-Nct1z 


<h  (z  ): 
P  P 


<n  (0)>  e 


=  <ii  (o)>  e 
P 


(3-30) 


The  right  hand  equality  is  a  statement  of  the  standard  definition 

of  the  scattering  cro  "s-section  of  a  material  with  N  scatters  per 
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unit  volume.  The  superscript  1  in  cr  indicates  that  this  is  the 

cross  section  for  a  single  mode,  and  integration  over  the  total 

number  of  modes  is  required  to  obtain  the  total  cross-section  (cr). 

-7  -] 

Since  Ncr  typically  has  a  value  of  about  10  cm  ,  <h^(z^)>  «  <hp(0) 

and  <n  >  can  be  treated  as  a  constant. 

P 

From  (3.3O),  we  get 


1 

a  = 


K(n.  +  1) 
N  c/Tlp 


and  comparison  of  this  to  (3.27)  yields 


1  Ss(n.  + 

c  - - 

N  <h  >  “n 
p  's 


(3.31) 


The  solution  for  n  (z  )  is 
s'  s' 


K(n.  +  l)<3i  > 
<h  (z  )>  =  - - - 

S  S  c/ti 


(3.32) 


Thus,  the  number  of  photons  (in  a  single  mode)  at  oo,  grows 

linearly  with  distance.  Using  representative  values  of  the  parameters 

-6 

in  (3.32)  leads  to  <h  >  =«  10  for  a  1  cm  laig  crystal.  This  justifies 

s 

the  assumption  <hg>  «  1  . 


Finally,  the  solution  for  n^(z^)  is 


Oi(zi)>  = 


,K<h  > 

_ £ _ 


a.  v.  -  K<n  > 
11  p 


a.  v.  -  K<n  > 
11  p 


1  -  e 


(3-33) 


For  spontaneous  processes,  v^  »  K  <hp>  j  this  reduces 
(3-33)  to 


<ft.(z.)>  = 

iv 


K<n  > 


a.  v. 
i  i 


-a.z.  \ 

i  i  • 


ii 

!  -  e  J 


(3-34) 


2  S  *1 

Also,  since  typically  =  10  -  ICr  cm  ,  (3-3^)  indicates 


that 


for 


K  <n  > 

<h.(z.)>  =  - =  constant 

1  1  a.v. 

i  i 


_2 

z.  >  ~  10  cm 

i 


(3-35) 


Thus,  the  number  of  photons  at  builds  up  to  the  value 

given  in  (3*35)  in  a  very  short  distance  and  remains  constant  at 
this  level.  Increasing  the  length  of  the  crystal  beyond  several 
"absorption  lengths"  (where  the  "absorption  length"  =  1/<X  )  does 
not  increase  the  number  of  idler  photons.  Physically,  this  result 
is  due  to  the  fact  that  at  the  value  of  <n^>  given  by  (3.35)  the 
absorption  and  generation  rates  of  idler  photons  are  equal  and  hence 
there  is  no  further  change  in  the  idler  photon  number. 

Substituting  representative  values  in  (3.35),  we  get 
^i(*i)>  ^  10  ^  ,  justifying  the  assumption  «  1  . 

The  results  of  this  section  are  plotted  in  Fig.  3.2. 
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FIG.  3-2--The  variation  of  photon  numbers  of  punp,  signal  and  idler, 
as  a  function  of  distance,  for  spontaneous  scattering.  The 
signal  photon  number  grows  linearly  with  distance,  while  the 
idler  photon  number  reaches  a  constant  value  after  a  short 
interaction  distance. 
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3*  Applications 


Using  the  results  of  the  previous  sections,  expressions  relating 
the  stimulated  gain  coefficient  and  spontaneous  scattering  parameters 
will  now  be  derived.  These  will  be  used  in  the  calculation  of  the 
material  nonlinear  coefficients  and  in  comparing  the  results  of 
spontaneous  scattering  measurements  to  the  developed  theory. 


a)  Relationship  of  Scattering  Cross-section  to  Stimulated  Gain 

The  total  cross-section  per  unit  solid  angle  and  per  unit 
frequency  range  (dcr/d^do)  is  obtained  by  writing 


do 

dfido 


=  cr 


±  dp(o>) 


dQdo 


(3-36) 


where 


dp(<u) 


dfldo 


is  the  number  of  modes  per  unit  solid  angle  per  unit  frequency  range, 
and  is  given  by 


dp(cu) 


2  2 
V  /(at  TC 


dn  do  (2it)^  '■ 


(3-37) 


Using  the  result  for  cr  (3.3I),  we  get  after  using  (3.27), 


dg  k(nt  -f  1)V  03  V2S  _  gs(Hj  +  1)  yig  03  V 
dfl  do  N(c/Tlp)(2ff)3  c3  Nc3<hp> 
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and  substituting  (3.25)  we  obtain  the  result 


da  g„  k"  hcu  (n  +  l) 

o  S _ p  1 _ 

dn  (2n)3  N  I 


(3.38) 


since 


This  is  the  derivation  of  the  result  quoted  as  (36)  in  HG. 
The  total  cross-section  is  obtained  by  integration  over  frequency 
and  solid  angle. 


b)  Calculation  of  S/Ldfid^ 

A  quantity  which  is  often  measured  is  the  ratio  of  scattered 
signal  power  to  pump  power  per  unit  crystal  length  per  unit  frequency 
range  and  per  unit  solid  angle;  that  is. 


P  Ldfido)  Ldfidu.' 
P 


where 

P 

S  =  —  .  (3.39) 

P 
P 


From  (3.3I)  and  (3*32),  we  get 

<n  >  Nti  da 

s  _ s  _ 

Tt  P  dT  da' 
P 


Ldnda.' 


Therefore, 


Ldn&a 


And  substituting  (3.38), 


O)  dcr 

—  N  - 

03  dfldo3 

p 


s 

gs  ft  cu3  (n.  +  1)1)2 

Ldncko 

(2*)3  ID  c2 

(3.40) 


This  shows  the  familiar  03g  dependence  of  the  scattered  power 


in  the  high-loss  region,  (i.e.  near  resonance)  since  in  this  region 
gg  a  <ds  [(3.23)  and  (3.27)]  . 

c)  Calculation  of  S/Ldfl 


A  convenient  experiment  is  the  measurement  of  right-angle 

scattering  (i.e.  9  -  90°).  The  results  of  these  measurements  are 

.  1*6 

given  in  terms  of  a  quantity  S/LdfJ  ,  where  integration  over  frequency 
of  S/Ldftcki?  is  performed  by  assuming  a  Lorentzian  lineshape  for  the 
mode  of  interest.  This  will  now  be  applied  to  (3. 40): 


S 

LdO 


/ 


MODE 


gs  b  oj  (r^  +  1)^ 

(2n  >3  Ip  c2 


(3-U) 


Making  the  assumption  of  a  "narrow"  mode  (that  is,  00  ,  T> 

s  s 

and  nj  are  assumed  constant  for  the  range  of  integration),  (l.4l) 
becomes 


S 

Ldf> 


h  cu3  ~2  (n.  ♦  1) 
(2 *)3  Ip  cc 


/ 

MODE 


*5 
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This  integrates  to 


JL  =  ft  4  ^(ni  +  l)r  gs(<D0) 

LdO  I6n2  I  c2 

P 


(3-42) 


since 

gs  =  Tg 

for  a  Lorentzian,  and)  £(<x>n  ,  m.  )  dm.  =  1.  ^ 

J  0  1 1  48 

Note  that  this  agrees  with  the  previously  quoted  result. 
(The  linewidth  there  is  defined  as  2 F  ;  hence  the  appearance  of 
the  numerical  factor  of  8  in  the  cited  reference. ) 


4.  Calculation  of  Spontaneously  Scattered  Idler  Power 

The  idler  power  (per  unit  pump  power,  solid  angle  and  frequency 
range)  generated  in  spontaneous  scattering  will  now  be  calculated. 
Using  the  expressions  for  idler  and  signal  photon  numbers,  (3.32) 
and  (3*35),  and  (3.25),  it  follows  immediately  that 


Fi  .  (aJ)  (— 

Pp  dn  dm  \0d  ms  /  V^dfidm 


A  certain  amount  of  care  must  be  exercised  in  the  application  of 
this  equation  to  an  experimental  situation.  This  is  discussed  in 
detail  in  Chapter  V  and  in  Appendix  E. 


(3-43) 


5*  Calculation  of  the  Nonlinear  Coefficients 


This  section  is  devoted  to  the  development  of  a  method  for  the 
calculation  of  the  nonlinear  coefficients  (d^  and  d^)  of  a  material. 
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The  method  involves  the  calculation  of  d^  and  d^  separately 
the  former  is  evaluated  from  right-angle  scattering  measurements; 
while  the  latter  is  computed  from  the  measured  second-harmonic 
generation  coefficient.  It  is  then  shown  how  these  computed  co¬ 
efficients  are  related  to  the  low-frequency  electro-optic 
coefficient. 


a)  Calculation  of  d„' 

- - - — tj 

To  compute  d'  ,  the  derived  expressions  for  the  on-resonance 


stimulated  gain  are  equated.  Starting  with  (3.15c) 

2  1  .  , 2  .  ,2 

m  0  I A  d' 

a  U  )  -  -5— £. — _ a. 

W  -  r 

V  “0  r 

and  substituting  (3.3)  applied  to  the  pump  wave 


1.  e. 


2n  I 


ell. 


we  get 


2tt  a)  S.  I„ 

«.K)  -  —  ’  0  f 

c  T)  7]  r 
s  ‘p 

2  2 

where  cjQq  has  been  used  (see  discussion  before(2. 12) ).  But 

from  (3. 1*2)  , 


gs(V  - 


I6n2  c2  (S/Ldp )  I 

b  ci>g  T|2  (n±  +  l)r 


therefore, 


A  *  + 

Snc1  T]  (S/Ldfl) 

1/2 

Q  “  “ 

_Sjto0“s  Vni  +  1}_ 

(3-44) 
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(Note  the  similarity  of  this  expression  to  (2.1 2), 


e  =  e  +  /  S . 

DC  o.  /.  J 


But.  from  (3* 9c), 


and  from  Butcher,' 


Low  freq 


"  2  ^Low 


*Low  freq  g^  ’ 

where  r  is  the  electro-optic  coefficient  and  T|  is  the  optical  re¬ 
fractive  index.  Therefore, 


d '  = 


"e  +  1  SJ  4<L  -  rT| 


(3.47) 


Thus,  the  calculated  values  of  and  d^  can  be  compared  to 

the  measured  electro-optic  coefficient  through  ( 3- ^7 ) • 

An  alternate  method  for  evaluating  the  ratio  of  cL'/d '  is 
.  37  «5i  r*  Q 

given  in  the  literature.  '  This  method,  which  compares  the 

intensities  of  spontaneous  scattering  at  the  transverse  resonant 

frequency  and  the  longitudinal  resonant  frequency,  was  first  applied 

to  GaP.  While  this  approach  works  well  for  a  material  with  a  single 

mode  (such  as  GaP),  the  mathematics  associated  with  applying  it 

to  a  multimode  material  is  rather  complicated,  and  the  solutions 

which  are  obtained  are  of  restricted  validity.  Therefore,  it  is 

not  used  in  this  work. 
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d)  Calculation  of  the  Raman  Susceptibilit; 


Using  the  previous  results,  the  Raman  susceptibility  is  calculated 
by  substituting  (3-44)  into  (3.9b);  this  yields 


c  TWs/Ldfl) 

V(»0)  =  — r^—z - 

0  2fio£  T)  (n.  +  l)r 


s  s  i 


(3-48) 


In  summary,  this  chapter  contains  the  analysis  of  the  stimulated 
and  spontaneous  scattering  processes.  The  first  part  is  devoted  to 
a  classical  analysis  of  stimulated  scattering.  The  important  result 
of  this  part  is  the  derivation  of  gain  expressions  for  multimode 
materials  which  are  valid  over  the  entire  mode.  The  latter,  half  of 


the  chapter  contains  a  quantum  mechanical  analysis.  The  spontaneous 
scattering  cross-section  is  calculated,  and  the  cross-section  is 
related  to  the  stimulated  gain  expression.  Methods  for  evaluating 
the  nonlinear  coefficients  are  developed,  and  an  expression  relating 
these  coefficients  to  the  electro-optic  coefficient  is  presented. 

We  now  turn  to  the  application  of  these  general  results  to  the 
specific  case  of  the  transverse  modes  in  LiNbO^,  since  the  experiments 
which  will  be  described  involved  interactions  with  these  modes. 


CHAPTER  IV 

APPLICATION  TO  THE  TRANSVERSE  MODES  IN  LiNbO, 

A.  INTRODUCTION 


In  this  chapter,  the  calculations  of  the  previous  chapter 
are  applied  to  the  transverse  modes  in  LiNbOy  Some  background 
material,  pertinent  to  this  work,  is  presented  briefly,  along  with 
the  dispersion  curve  for  the  A^-symmetry  modes.  Then  the  nonlinear 
coefficients  are  calculated  and  compared  to  the  measured  electro- 
opt.ic  coefficier  Finally,  the  scattering  parameters  for  the  two 
important  transverse  modes— the  248  cm  1  and  628  cm  ^  modes-  are 
caculated. 


B.  BACKGROUND 


LiNbO,  has  been  studied  extensively  and  many  of  its  properties 
■2  31  34  52  53 

are  described  in  the  literature.  ’  A  brief  summary  of  the 

properties  important  for  this  work  is  given  below. 

LiNbOj  is  a  uniaxial  crystal  in  point  group  C  ,  with  the 

axis  of  symmetry  defined  as  the  c-axis.  If  a  light  beam  traveling 

through  the  material  is  polarized  perpendicular  to  the  c-axis  (E  ]_  c), 

the  appropriate  index  of  refraction  is  t]q  ,  the  ordinary  index.  If 

the  beam  is  polarized  so  that  E  ||  c  ,  -rj  =  -rje  ,  the  extraordinary 

index.  Similarly,  the  normal  modes  of  vibration  of  the  crystal  are 

classified  according  to  their  symmetry  with  respect  to  the  axes  of 

symmetry  (the  c-axis  and  the  three  mirror  planes).  Vibrations  which 

are  symmetric  with  respect  to  both  the  c-axis  and  the  mirror  planes 

are  classified  as  A  -modes.  The  existence  of  infrared  and  Raman  activity 
1  54 

can  be  determined  from  Hertzberg  and  it  is  found  that  the  A^  modes 
are  both  infrared  and  Raman  active.  This  ensures  the  existence  of 
dj/  and  d^'  and  the  possibility  of  producing  tunable  signals  in 
interactions  with  these  modes. 
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The  appropriate  dispersion  curve  is  the  one  obtained  for  E  jj  c  , 
since  the  properties  of  this  curve  are  determined  by  the  A^-symi  try 
modes.  All  of  the  quantities  needed  to  evaluate  the  propagation 
characteristics  have  teen  measured  by  Barker  and  Loudon^  and  are 

tabula' ed  in  Fig.  4.1. 

Plots  of  k'.  and  k'.'  as  functions  of  v.  were  obtained  from 
11  1 

Eq.  (2.1l),  where  and  k  ( "  are  the  real  and  imaginary  parts  of 

the  propagation  constant  respectively.  These  computer-genera4' °i  curves 

are  shown  in  Fig.  4.2.  The  two  large  "bumps"  on  the  kC  plot  are 

-1  - 1  "L 

associated  with  the  248  cm  and  628  cm  vibrational  frequencies. 

Therefore,  the  region  for  y  .  «  2^ 8  cm  ^  is  called  the  "248  cm  ^  mode" 

_2  _ 

and  similarly  for  v  ^  551  628  cm  ,  the  region  is  called  the  "628  cm 

mode."  Consideration  of  the  scattering  from  these  two  modes  is  the 

subject  of  the  remainder  of  this  study. 

The  wave  vector  matching  condition  is  now  imposed.  The  solution 

of  Eq.  (2.l4)  (see  Fig.  2.4)  yields  a  family  of  phase-matching  lines. 

(Note  that  these  lines  depend  upon  the  choice  of  pump  frequency.) 

The  values  for  r\  and  q  were  taken  from  the  data  of  Hobden  and 
55  P  s 

Warner,  using  the  equation  for  qg  .  The  computer  generated 
phase-matching  lines  for  a  ruby  laser  pump  (\  =  6945  X),  super¬ 

imposed  on  the  -symmetry  dispersion  curve,  are  shown  in  Fig.  4. 5, 

C.  CALCULATION  OF  THE  NONLINEAR  COEFFICIENTS 

In  this  section,  the  calculations  of  the  nonlinear  coefficients  for 
the  A^- symmetry  modes  of  LiNbO^  are  presented.  The  results  are 
restricted  to  interactions  in  which  all  waves  are  polarized  parallel 
to  the  c-axis,  (extraordinary  waves) .  This  conforms  to  the  experimental 
situation  described  in  Chapter  V. 

Therefore,  the  pertinent  nonlinear  constants  are  those  with 

*  .  .  56 

subscript  j>33  (or  33  in  contracted  notation). 

In  the  convent. ional  notation,  the  z-direction  is  taken  along  the 
o-axis,  and  this  direction  is  also  denoted  by  the  numeral  3*  Hene»>, 
d-. .  is  the  appropriate  constant  when  all  fields  and  polarizations  ar-* 
at 6fig  the  c-axis. 
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FIG.  4.1— Linear  properties  of  A^-synmetry  node'  of  LiNbO,,  as 
determined  by  infrared  reflectivity  measurements  of 
Barker  and  Loudon.  lhe  mode  at  692  cm‘^  has  been 
identified  as  a  combination  band,  since  group  theory 
predicts  the  existence  of  only  four  fundamental 
symmetry  modes.  Its  contribution  to  the  dielectric 
constant  must  be  included. 
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FIG.  <*.2--The  right  side  of  the  figure  is  the  dispersion  curve  for  the  A^-syjssnetr;,- 
modfcs  of  LiirbO^;  the  left  side  shews  frequency  vs  electric  field  loss.  The 
curves  were  Renerated  from  a  computer  progress,  using  the  dnta  of  Fig.  It.1  . 
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From  Eq.  (3.45), 


E 


d3  =  i6Hi._  =  5 .28  <  io  ■  -‘p .  , 

jj 


where  the  latest  published  measurement  of  d  has  beer,  us- 


Q  ' 


From  Eq.  (3.44), 


dQ' 


tjffC  np(s33/^a) 


t/  ■  _ - 

S  .ftuo  o>  n  (n.+l) 
j  0  s  >s  l 


1/2 


esu 


The  \alues  of  S  /lid£  for  the  modes  of  LitfbCL  i.avn  Iv¬ 
lin  ,  3 

measured  and  the  resulting  values  for  <dA  are  tabulated 

#  W 

Fig.  4.4. 


Raman  susceptibility. 

The  Raman  sv;sceptibility  is  calculated  from  Eq .  (3.4.?,)  vein.,- 
the  measured  values  of  S  /  I4H  .  The  results  are  tabulate!  ir. 
Fig.  4.1*, 

Finally,  there  results  are  compared  to  the  measured  ejeetro- 
optic  coefficient  through  Eq.  (3-47).  Turner  has  rrr<rte5'"J 

"  3°.8  x  10‘10  cm/v  -  *  !0'(  -.-u  , 


There  is  ar.  arbitrariness  ir  the  determination  r  th-  r  t*r. 
dq  .  The  approach  of  Kami  now  and  Johnston  °  will  be  f«. . 
regard. j  This  leads  to  negative.-  values  of  dq  for  the  ?"*4  .«*’*  a 
30?  cm  modes.  It  should  be  pointed  out  that  Gap  has  a  r-gv  -  - 
value  of  d^  lor  equivalently  •  ratio  of  (d£  /do)  <  0  ’.>*  tr. 
addition,  it  should  also  be  noted  that  arbitrarily  chore  1  tv  *.-«»  r 
of  each  <iA  to  be  positive  would  increase  the  calculated  v*ootr< 
T*  ic  coefficient  by  only  about  . 
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4. 3- -The  phase-matching  lines  for  a  ruby  laser  pump  (x=  69^3  &)  are  super¬ 
imposed  on  the  dispersion  tCurve  of  Fig.  b.Cl.  The  optical  refractive  index 
data  of  Hobden  and  Warner  ?  were  used . 
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Fig.  k.k — Tabulated  properties  of  the  Aj_  -symmetry  modes  of  LiNbO  .  The  values  of  <i'  a  '  and 
XR  (v  )  are  the  significant  results  of  this  section.  ^  ^ 


and  therefore 


4  -5 

tj  r^  =  2.36  x  10  esu 

A  comparison  of  this  result  to  the  calculated  value  of 
h",  +  Z  S-d,/  (Fig.  4.4)  shows  the  two  to  be  in  excellent 

agreement.  J 

D.  SCATTERING  PARAMETERS  FOR  THE  248  CM-1  AND  6 28  CM-1  MODES 

Using  the  constants  calculated  in  the  previous  section  and  the 
development  of  Chapter  III,  we  devote  this  section  to  the  computa¬ 
tion  of  the  scattering  parameters  for  the  248  cm  ^  and  628  cm  1 
transverse  modes  in  LiNbO^.  With  reference  to  Fig.  4.2,  several 
comments  are  in  order  before  proceeding  with  the  calculation.  These 
comments  are  intended  to  point  out  two  distinctions  between  the  modes, 
which  are  pertinent  to  these  calculations: 

a)  The  absorption  (a.  =  2  k.")  along  the  628  cm-1  mode  is 

f  1  2  -1 

high  over  the  entire  mode,  varying  from  l4  x  10  cm  at  the 

o  2  1  o 

0=0  point  to  200  x  10  cm  for  phase-matching  at  6  .  Since 

stimulated  gains  are  of  the  order  of  1-10  cm  1  for  pump  powers 

available  from  standard  Q-switched  lasers,  the  idler  loss  is  always 

much  greater  than  the  gain  and  the  high  loss  approximation,  Eqs. 

(3.15)  cr  (3.18b),  is  valid  over  the  entire  mode. 

The  high-loss  approximation  does  not  apply  to  the  case  of  the 

243  cm  ^  mode  however,  since  the  loss  varies  from  approximately  zero 
2  -1 

to  50  x  10  cm  over  the  near- forward  scattering  region.  Hence, 
the  more  general  expression  is  required. 

b)  Referring  to  the  9=0  line  in  Fig.  4.3,  it  is  seen  that 

collinear  scattering  from  the  628  cm  mode  is  possible  and  indeed 
22  23  -1 

has  been  observed.  *  For  the  248  cm  mode,  on  the  other  hand, 
collinear  forward  scattering  is  not  obtainable,  since  the  intersection 
of  the  6=0  line  occurs  at  V  =  0  .  Thus,  the  collinear  gain 
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expression, Eq.  (J.l6),  is  ncc  strictly  applicable  at  any  point  on 
this  mode. 

Using  calculated  coefficients,  numerical  values  for  the  signal 

stimulated  gain  (g  )  were  computed  for  the  two  modes  and  the 
s 

results  are  plotted  in  Fig.  4.5.  These  curves  are  for  a  pump  power 
8  2 

density  of  5  X  10  W/cm  ,  typical  of  focused  Q-s witched  lasers  and 
the  approximate  power  density  used  in  the  experiments  to  be  described 
later.  Note,  that  for  all  points  on  the  dispersion  curve,  except  those 
very  near  to  V.  =  0  ,  g  is  greater  for  scattering  from  the 

-1  ^  ®  _ i 

248  cm  mode  than  for  the  628  cm  mode. 

Figure  4.6  shows  the  calculated  idler  power  (divided  by  signal 
power)  as  a  function  of  frequency  for  stimulated  scattering  from  the 
two  modes.  Again,  it  is  noted  that  in  scattering  from  the  243  cm-1 
mode,  the  ratio  P./P  exceeds  that  from  the  628  cm-1  mode  by  more 
than  an  order  of  magnitude. 

Similar  plots  for  the  spontaneous  case  are  given  in  Figs,  4.7 

and  4.8  where  P  /P  L  dfldo)  and  P./P  dfldcD  are  plotted  for  both  modes, 
s  p  ^-P 

As  for  the  stimulated  process,  the  scattered  power  (both  signal  and 
idler)  is  greater  for  the  248  cm  mode  than  the  628  cm"1  mode.  The 
shape  of  the  curve  for  the  628  cm-1  mode  in  Fig.  4.7  has  been  verified 

17 

expermentally  by  Puthoff,  et  al,  whose  measurements  indicated  a 
nearly  constant  signal  power  with  frequency  in  spontaneous  scattering 
for  this  mode. 

In  the  next  chapter,  we  shall  compare  the  theoretical  predictions 
to  the  experimental  results.  The  agreement  between  theory  and  experi¬ 
ment  will  be  seen  to  be  good  for  both  stimulated  and  spontaneous 
scattering . 
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FIG.  4.6--Idler  power  (divided  by  signal  power)  generated  in  stimulated 

scattering  frop  the  2U8  cm"1  mode  (Curve  a)  and  the  628  cm-1  mode 
(Curve  b).  Note  the  use  of  different  ordinates  for  the  two 
modes.  Equation  (3.26)  and  the  calculated  constants  of  Fig. 

4.U  were  used  to  plot  the  curves. 


ropup  d 


FIG.  U.8 — Idler  power  (per  unit  pump  p  vor,  o  lid  angle  end  bandwidth) 

fenerated  in  spontane-UG  scattering  f'r  ir.  the  ■M  cm"1  rav  dr 
Curve  a)  and  the  628  mb"1  m  ie  (Curve  l<).  The  rdinato  •« 
in  logarithmic.  For  oroail  v.  ,  Curve  a  "turns  ar- und"  ats  i 
intersects  the  origin;  thi6  nas  been  omitted  for  clarity. 
Equation  (3.**3)  and  the  calculated  constanto  of  Fig.  ** . *»  were 
used  to  plot  the  curveG. 


CHAPTER  V 


EXPERIMENTAL  EFFORT 


A.  INTRODUCTION 

The  experimental  effort  was  divided  into  three  areas: 

a)  The  first  set  of  experiments  was  an  attempt  to  detect  the 
idler  ( infrared)  power  generated  in  spontaneous  scattering  from  the 
628  cm  mode.  An  argon- ion  laser  was  built,  and  a  phase-sensitive 
detection  system  was  assembled.  This  experimental  set-up  is 
adaptable  to  a  wide  variety  of  scattering  measurements.  The  experiment 
was  not  definitively  successful,  in  that  the  detection  of  the  idler 
could  not  be  conclusively  verified. 

b)  Next,  a  series  of  stimulated  scattering  experiments  using  a 
Q-s witched  ruby  laser  were  performed.  The  most  significant  experiments 
yielded  tunable  stimulated  optical  emission  in  scattering  from  the 

248  cm  1  mode.  No  external  resonator  was  required  to  produce  the 
tunable  output.  In  addition,  the  idler  wave  involved  in  this  process 
was  detected.  All  of  the  measured  parameters  agree  reasonably  well 
with  the  theory  presented  in  the  previous  chapters. 

c)  Finally,  measurements  of  spontaneous  scattering  from  the 
248  cm  1  mode  were  made.  Both  the  frequency  and  magnitude  of  the 
scattered  signal  as  a  function  of  phase-matching  angle  6  were 
measured  and  found  to  agree  well  with  prediction. 

B.  INFRARED  GENERATION  BY  SPONTANEOUS  SCATTERING  FROM  THE 

628  cm’1  MODE. 

At  the  start  of  the  experimental  program,  tunable  stimulated 
Raman  scattering  had  not  as  yet  been  observed.  Measurements  of  the 
tunable  signal  in  spontaneous  scattering  from  the  366  cm’1  mode  in 
GaP, the  Ul4  cm’1  mode  in  ZnO,^  and  the  628  cm’1  mode  in  LiNbO^,1^ 
had  recently  been  performed.  Starting  with  the  results  of  the  last 
work,  the  objective  of  the  first  part  of  the  experimental  program 
was  the  detection  of  the  tunable  idler  power  generated  in  spontaneous 
scattering  from  the  628  cm’1  mode  in  LiltbO^.  The  two  quantities 
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of  interest  are  the  frequency  tuning  and  the  idler  power. 

The  tuning  curve  is  obtained  from  the  computer  generated  data  used 
to  plot  Fig.  4.3*  Since  the  three  wave  vectors  are  available  for  each 
value  of  6  ,  the  angle  <p  is  calculable.  The  result  is  plotted  in 

Fig.  5.1.  Note  that  cp  is  large  over  most  of  the  tuning  range.  Due 
to  the  large  index  of  refraction  at  the  idler,  most  of  the  infrared 
power  would  undergo  total  internal  reflection  upon  striking  a  surface 
normal  to  the  incident  laser  beam.  For  this  reason,  the  crystal  was 
cut  at  ^5° at  the  output  face  as  shown  in  Fig.  5*2.  The  calculated  external 
tuning  curve  is  also  shown  in  this  figure  and  it  is  seen  that  a  tuning 
range  of  525  cm  1  to  595  cm  1  (17  p  to  19  p)  is  expected. 

The  scattered  infrared  power  was  calculated  from  the  results  of 
Chapter  IV,  with  particular  reference  to  Fig.  ^.8,  where  P^/P  dfjdm 
is  plotted  as  a  function  of  .  However,  we  are  interested  in  the 

scattered  idler  power  in  a  specified  frequency  interval  and  solid 
ang3°  as  a  function  of  cp"  (defined  in  Fig.  5*2).  This  is  shown  in 
Appendix  E  to  be: 


In  the  above  expression: 

is  the  idler  power  in  a  frequency  range  scattered 

into  an  aperture  of  height  H  at  a  distance  r  from  the 
crystal  output,  £$'  is  the  angle  increment  corresponding 
to  Ca>  (see  Fig.  5*2)  , 

and 

d<p/de  is  calculated  from  Fig.  5*1* 

A  sketch  of  the  geometry  is  given  in  Appendix  E.  Carrying  through 

11 

this  calculation  for  —  »  0.1  and  P  ■  0.J  H  we  get  a  total 

r  P  -12 

power  over  the  tuning  range  of  about  3  x  10  watts.  If  the  range  is 
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l/\  (cm 


^  (cm1) 


divided  into  six  segments,  (each  encompassing  a  &p*  of  about  15 
and  a  of  10  cm  ),  the  power  in  each  segment  is  about 

0.5  x  lO"-^  watts.  As  will  be  described  later,  this  is  only  slightly 
above  the  limit  of  detectability  of  the  phase  sensitive  detection 
system  employing  the  Ge:Cu  detector.  Even  though  the  probability 
of  success  of  the  experiment  appeared  marginal,  the  significance 
of  the  result  (if  %hievable)  justified  the  experimental  attempt. 

In  addition,  the  calculated  tuning  curve  is  valid  for  stimulated 
scattering,  and  would  be  useful  in  such  experiments. 

IT  52  59 

Most  of  the  previously  reported  scattering  measurements  *  7 
used  a  He-Ne  laser  (6328  X)  as  the  pump  source.  The  use  of  an  argon 
laser  instead  has  the  following  advantages: 

a)  The  argon  laser  has  a  power  output  cf  about  an  order 
of  magnitude  greater  than  the  highest  power  conventional  He-Ne 
laser.  It  is  comparatively  easy  to  get  an  output  of  520  mW  at 
4880  X  using  the  argon  laser,  while  100  mW  is  the  output  of  the 
highest  power  commercially  available  He-Ne  laser.  Since  the  scattered 
power  is  proportional  to  pump  power,  this  advantage  is  obviously 
significant. 

b)  The  shorter  wavelength  is  also  an  advantage  for  two  reasons. 
First,  scattered  power  is  proportional  to  the  fourth  power  of  a 
signal  frequency  [see  (5.4o)].  This  yields  about  a  factor  of  3 
increase  in  scattering  efficiency  for  a  4880  X  pump  compared  to 

a  6328  X  source.  Secondly,  since  most  phototubes  are  more  sensitive 
in  the  green-blue  spectral  region  than  in  the  red,  it  is  possible 
to  detect  weaker  scattered  optical  signals  with  a  4880  X  pump. 

c)  Finally,  several  wavelengths  are  available  from  an  argon 
laser  ( including  the  two  strong  lines  at  4880  X  and  5145  X) .  This 
mny  be  useful  in  certain  experiments . 

For  these  reasons,  the  argon  laser  was  selected  as  a  pump  source 

We  also  decided  to  construct  the  l'ser  in-house,  employing  rf 
60 

excitation,  and  incorporating  several  features  not  available  in 
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commercial  units: 


(1)  A  self  contained  vac  mm  system,  including  both  roughing 
and  Vac-ion  pimps.  This  system  allows  for  adjustment  and  monitoring 
of  operating  gas  pressure  as  well  as  regular  "flushing"  of  the  system. 

(2)  Isolated  supplies  of  argon  gas.  This  permits  operation 
of  the  system  with  other  gases,  e.g.,  krypton. 

j  ( 3)  Ball-socket  construction  of  end  windows  allows  for 

i  relatively  easy  cleaning  and/or  replacement  of  Brewster  angle  windows. 

(4)  Interchangeable  parts.  A  spare  laser  tube  was  built  and 
|  was  easily  installed  into  the  system  after  failure  of  the  initial  tube, 

j  ^  5)  Mirror  flexibility.  The  external  mirror  mounts  provide 

1 

!  simplicity  of  alignment  and  adjustments  of  the  optical  cavity,  as  wall 

as  easy  interchangeability.  The  system  operates  with  a  6m  spherical 
{  ~5$  transmission)  output  mirror.  The  back  mirror  is  either  a 
| .  totally  reflecting  flat  which  yields  output  at  all  wavelengths 

simultaneously  or  a  prism  with  a  totally  reflecting  back  face  which 
allows  lasing  action  at  a  single  wavelength  only.  In  this  latter 
i  configuration,  the  output  wavelength  can  be  changed  from  one  of  the 

lasing  lines  to  another  simply  by  rotating  the  prism. 

In  addition,  the  in-house  construction  of  the  laser  resulted  in 
a  working  system  at  considerably  less  expense  than  comparable 
commercially  available  lasers.  A  photograph  of  the  laser  system  is 
shown  in  Fig.  5. 3*  (An  aluminum  box  was  built  around  the  laser, 

subsequent  to  the  taking  of  the  photograph,  to  reduce  pickup  in 
the  detection  system.) 

The  laser  performance  was  very  close  to  the  design  goals,  with 
a  total  power  output  of  about  1.25  watts;  the  power  was  divided  among 
the  various  wavelengths  as  shown: 


Wavelength  ($) 

Power  (MW) 

4727 

50 

4765 

130 

4880 

600 

4965 

100 

5018 

20 

5145 

340 

FIG.  5. 3- -Photograph  of  the  argon  laser  system.  The  excitation  is  supplied 
by  a  3  MHz  oscillator  (shown  in  background  center).  The  interchangeable 
Brewster's  angle  windows  are  at  the  ends  of  the  tube.  The  coated  prism  is 
in  the  mirror  mount  at  the  far  left  of  the  table;  the  beam  emerges  from  the 
output  mirror,  which  is  at  the  far  right  of  the  table.  The  system  is  evacu¬ 
ated  with  a  roughing  pump  (not  visible)  and  a  Vac-ion  pump  (near  center  of 
photograph) .  The  argon  gas  is  stored  in  the  flasks  seen  in  the  foreground; 
the  hoses  are  for  water  cooling  of  the  tube. 


The  limitation  on  the  power  output  of  the  laser  was  imposed  by  the 
oscillator:  its  maximum  output  was  about  5  KW. 

Detection  of  the  scattered  signals  was  accomplished  using  a 
phase-sensitive  detection  system,  employing  a  lock-in  amplifier 
(Princeton  Applied  Research  Corp.,  Model  HR-8).  The  main  features 
of  this  technique,  which  is  described  in  detail  in  the  literature, 
are  summarized  below: 

The  signal,  which  is  chopped  at  a  fixed  frequency,  is  amplified 
in  a  low-noise  signal  amplifier  whose  center  frequency  and  bandwidth 
are  variable.  Simultaneously,  the  chopper  generates  a  reference  signal 
(by  means  of  a  light  bulb  and  phototube)  which  is  amplified  in  a 
reference  channel.  The  reference  and  signal  are  then  mixed  and  a  dc 
output,  proportional  to  the  magnitude  of  the  signal  and  the  cosine 
of  the  phase  angle  between  the  signal  and  reference,  is  obtained. 

This  is  th  .1  passed  through  RC  low  pass  filtering. 

The  primary  advantage  of  this  system  is  that  it  can  extract  weak 
signals  from  the  noise  background.  This  is  accomplished  primarily  by 
narrowing  the  bandwidth  of  the  signal  channel  to  eliminate  white  noise 
and  by  synchronously  comparing  the  signal  and  reference. 

The  detector  employed  in  the  experiment  varied  with  the  type  of 
measurement  being  made.  For  optical  measurement,  a  photomultiplier 
tube  was  used,  while  for  the  infrared  measurements,  a  Ge:Cu  semi¬ 
conductor  detector  cooled  to  liquid  helium  tempera  re  was  employed. 
This  is  described  in  more  detail  below. 

The  experimental  apparatus  is  shown  in  F^g.  5.4.  The  laser  beam 
was  chopped  at  a  1  kHz  rate,  while  a  reference  signal  was  simultaneously 
generated.  The  laser  beam  then  entered  the  LiNbO  crystal  as  an 
extraordinary  wave  and  was  totally  internally  reflected  at  the  45  cut 
in  the  output  face.  The  acceptance  angle  of  the  idler  output  was 
determined  by  an  aperture,  and  the  detected  band  of  idler  frequencies 
was  varied  by  translation  and  rotation  of  the  adjustable  mirror.  After 
passing  through  a  J  meter  Jarrel  Ash  monochromater,  which  had  a  grating 
blazed  at  18  ji,  the  signal  was  detected  in  the  Ge:CU  detector;  the 
output  of  this  detector  was  fed  into  the  signal  channel  of  the  lock-in 
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FIG.  5 .4- -Experimental  apparatus  for  detecting  scattered  idler  radiation. 
The  laser  (shown  schematically  as  a  tube,  T  ,  a  mirror,  ,  and  a 

coated  prism,  P)  output  beam  (about  500  raW  of  power  at  4880  A)  was  chop¬ 
ped  at  a  1  kHz  rate  by  the  chopper,  C  ,  while  a  reference  signal  was 
generated  simultaneously  as  the  photodetector,  PD  ,  detected  the  chop¬ 
ped  output  of  a  light  bulb,  LB  .  The  beam  passed  through  the  LiNbOj 
crystal,  XTAL  ,  as  an  extraordinary  wave  (the  laser  polarization  and 
the  c-axis  of  the  crystal  were  both  normal  to  the  paper) .  The  laser  beam 
was  totally  internally  reflected  at  the  45°  cut  in  the  output  face.  The 
idler  radiation  was  reflected  by  the  rotatable,  translatable  mirror,  , 
through  the  aperture,  A  .  It  then  passed  through  the  monochromator, 

MC  ,  and  was  detected  by  the  Ge:Cu  detector,  D  ,  which  was  cooled  to 
4° K  in  a  dewar.  The  detector  output  and  the  reference  signal  were  fed 
into  the  PAR  Model  HR-8  synchronous  detector,  SD  ;  the  output  of  the 
detection  system  was  monitored  by  a  chart  recorder,  R  . 
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amplifier.  Finally,  the  dc  output  of  the  amplifier  was  recorded  on 
a  chart  recorder. 

The  performance  of  the  detector  was  one  of  the  crucial  parameters 

affecting  the  success  of  the  experiemnt.  A  Ge:Cu  detector,  purchased 

from  Santa  Barbara  Research  Center,  a  subsidiary  of  Hughes  Aircraft 

Company,  was  used.  The  properties  of  this  type  of  detector  have  been 

discussed  in  the  literature. ^*^5  m0st  significant  features,  as 

far  as  this  experiment  was  concerned,  were  its  high  detectivity  and 

responsivity,  low  noise  equivalent  power  and  speed  of  response.  Its 

major  disadvantage,  as  compared  to  a  Golay  cell,  for  example,  is  that 

cooling  to  4.2°K  is  required.  Based  upon  the  detector  parameters 

supplied  by  the  manufacturer,  it  was  estimated  that  by  using  the  phase- 

sensitive  detection  system  it  would  be  possible  to  detect  a  signal  of 
-1^ 

approximately  10  watts.  This  is  only  a  factor  of  two  larger  than 
the  expected  power  output  through  the  aperture  in  Fig.  5.4. 

With  the  experimental  set-up  shown,  no  idler  signal  was  observed. 
Since  the  calculations  indicated  that  the  signal  power  was  marginal 
for  detection,  this  result,  although  disappointing,  was  not  too 
surprising.  Attempts  to  focus  the  signal  using  Ge  lenses  failed, 
(probably  because  of  line-up  inaccuracies).  Finally,  the  mirror, 
aperture,  and  the  monochromator  were  removed  and  the  crystal  was 
mounted  inside  the  dewar  at  a  distance  of  about  -J-  inch  from  the  detector. 
For  this  configuration,  the  solid  angle  viewed  by  the  detector  (and 
hence  the  signal  power)  was  increased  by  about  two  orders  of  magnitude, 
and  the  losses  in  the  path  from  the  crystal  to  the  detector  were 
eliminated.  Signals  were  now  observed  at  about  the  predicted  level. 
However,  it  was  then  ascertained  that  the  Ge:Cu  detector  was  sensitive 
to  visible  radiation.  Since  the  infrared  signal  wa^  predicted  to  be 
about  18  orders  of  magnitude  below  the  pianp  power,  it  was  difficult 
to  rule  out  the  possibility  that  the  detector  was  responding  to  stray 
pump  light,  rather  than  to  infrared  radiation.  Various  attempts  to 
prevent  the  visible  radiation  from  reaching  the  detector  were  then  made, 
but  the  results  were  inconclusive. 

In  summary,  it  could  not  be  conclusively  established  that  the 
detection  of  the  infrared  power  had  been  achieved. 
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C .  TUNABLE  OPTICAL  AND  INFRARED  GENERATION  BY  STIMULATED  SCATTERING 

FROM  THE  248  cm-1  MODE. 

22  23 

The  reports""  ’  of  stimulated  scattering  from  LiNbO^  led  to  a 

redirection  of  the  experimental  effort  to  the  study  of  tunable  infrared 

power  in  stimulated  scattering  from  the  628  cm-1  mode.  Early  in  the 

course  of  the  stimulated  experiments,  a  new  process -tunable,  stimulated 

optical  scattering  from  the  248  cm  ^  mode,  without  an  external  resonator- 

was  observed.  In  addition,  the  tunable  infrared  wave  generated  in  this 

24 

process  was  detected.  The  experiments  have  been  reported,  and  a 
more  detailed  account  has  been  given  by  Yarborough.^  The  highlights 
of  the  experiments,  particularly  those  features  which  relate  to  the 
theory  developed  in  this  work,  are  summarized  below. 

The  experimental  apparatus  is  shown  in  Fig.  5.5.  The  pump  was  a 
Q-switched  ruby  laser  emitting  ~  20  nsec  pulses  with  a  peak  power  of 
approximately  1  MW  and  a  spot  size  at  the  laser  of  about  2  mm.  A 
lens  focused  the  beam  into  a  3*3  cm  a-axis  crystal,  with  the  laser 
polarized  along  the  c-axis.  The  end  faces  of  the  crystal  were  polished 
flat  and  parallel  to  within  a  few  seconds  of  arc.  Temporal  behavior 
of  the  input  pump,  output  pump,  and  signal  was  recorded  with  photodiodes 
und  the  wavelengths  were  recorded  on  film  at  the  output  of  aim 
grating  spectrometer.  The  spatial  separation  of  the  beams  was  detected 
by  photographing  the  beams  on  a  screen. 

Figure  5.6  shows  the  output  wavelengths  and  spatial  distributions 
of  the  signal  for  two  different  angles  of  incidence,  4>  ,  between  the 

pump  beam  and  surface  normal.  A  50-cm  focal-length  lens  was  used  in 
both  cases.  With  this  focal-length  lens,  crystal  damage  was  not 
observed.  Figure  5.6(a)  shows  the  spectrometer  output  for  the  case 
of  normal  incidence  $  ■  0°  ,  and  5.6(b)  shows  the  corresponding 

spatial  separation  at  a  distance  of  5*6  m  from  the  crystal.  This 
separation  implies  an  angle  of  0,42°  between  the  signal  and  pump 
beams  ii.side  the  crystal.  The  observed  frequency  shift  in  this  case 
is  50  cm  ,  and  the  observed  spatial  separation  agrees  with  this  shift. 
This  point  is  plotted  as  point  A  on  the  dispersion  curve  in  Fig.  5.7. 
Figures  5.6(c)  and  5.6(d)  show  the  spectral  and  spatial  outputs  from 
the  crystal  for  an  angle  of  incidence  $  ■  2.30°.  The  observed 
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FIG.  5- 6- -Photographs  of  signal  wavelength  and  spatial  shifts.  The  wave¬ 
lengths  were  photographed  at  the  spectrometer  output  (a  and  c); 
the  spatial  positions  of  the  beams  were  photographed  at  a  dis¬ 
tance^  5-6  meters  from  the  crystal  (b  and  d).  a  and  b  are  for 
0=0  (normal  incidence);  c  and  d '’are  for  0  =>  2.3°.~  The 
wavelength  and  spatial  photographs  are  in  good  tgreement  with  the 
predictions  from  the  dispersion  curve  (see  Fig.  5.5).  For  these 
pictures,  the  laser  power  was  abouc  600  kW  and  was  focused  by 
a  50  cm  lens. 
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frequency  shift  in  this  case  is  13^  cm  \  and  again  the  spatial 
separation  agrees  with  the  observed  shift.  This  point  has  been 
plotted  as  B  on  the  dispersion  curve.  In  both  Figs.  5.6(b)  and 
5.6(d),  the  lower  beam  is  the  laser  and  the  upper  beam  is  the  signal. 
Other  measured  points  have  also  been  plotted  in  Fig.  5.7. 

Figure  5.8  shows  the  typical  temporal  response.  Figure  5.8(a) 
shows  the  incident  pump  pulse,  5.8(b)  the  transmitted  punp  pulse, 
and  5.8(c)  the  corresponding  signal  pulse.  For  this  particular  case, 
the  conversion  was  about  60$.  The  maximun  conversion  observed  was 
approximately  70$. 

In  Fig.  5.9  the  observed  signal  frequency  shift,  (and  corresponding 
idler  wavelength),  is  shown  as  a  function  of  the  angle  of  incidence 
of  the  laser  beam  relative  to  the  normal  to  the  crystal  surface.  The 
presence  of  the  parallel  crystal  faces  makes  possible  the  tuning  of 
the  signal  as  the  angle  of  incidence  is  varied.  When  the  crystal  was 
repolished  so  that  there  was  an  angle  of  1.5°  between  opposite  crystal 
surfaces,  it  was  no  longer  possible  to  tune  the  signal  by  varying  $  . 

In  addition,  it  was  necessary  to  increase  the  pump  intensity  in  order 
to  obtain  a  comparable  signal  intensity. 

Good  agreement  with  the  tuning  curve  of  Fig.  5*9  has  been 
obtained  by  the  "two-bounce"  theory  of  Yarborough.^  This  theory 
assumes  that  the  signal  beam  makes  two  forward  passes  through  the 
puap  beam  inside  the  crystal.  The  second  forward  pass  is  due  to  a 
reflection  from  the  output  crystal  surface  and  a  re-reflection  at.  the 
input  surface.  This  explains  the  role  of  the  parallel  crystal  faces 
in  the  tuning  of  the  signal  frequency.  By  using  the  appropriate 
crystal  dimensions  and  pump  beam  spot  site,  a  good  fit  is  obtained  tc 
the  measured  tuning  curve. 

An  estimate  of  the  gain  constant  can  be  made  from  the  experimental 

66 

parameters.  Tht  result  it  g  «  >  cm*  ,  as  swing  normal  incidence 

a 

of  the  pump  and  assuaing  that  the  "tvo-bource"  theory  is  valid. 

Direct  detection  of  the  infrared  idler  was  accomplished  with  the 
experimental  set-up  shown  in  Fig.  5.10.  As  shown,  a  60°  cut  was  made 
in  the  corner  cf  the  output  end  of  the  crystal  to  allow  the  idler  wive 
to  emerge  approximately  normal  to  the  crystal  surface,  thus  minimising 
reflection  loss  at  the  crystal-air  interface.  To  filter  out  stray 
optical  signals,  a  strip  of  blacli  polyethylene  was  placed  in  the  light 
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pipe  leading  to  the  InSb  detector,  cooled  to  1.5°K.  The  output 
of  the  detector  was  monitored  on  a  Tektronix  Ttype  555  oscilloscope. 
Figure  5.11  shows  the  idler  pulse  accompanying  a  pump-to-signal 
frequency  shift  of  ~5°  cm”1  ,  which  corresponds  to  an  idler  wavelength 
of  200  u  .  Based  upon  measurements  of  the  pump  and  signal  pulses 
using  photodiodes  [Figs.  5.8(a)  -5.8(c)],  the  idler  pulse  duration  is 
expected  to  be  about  20  nsec .  The  idler  pulse  stretching  may  be 
attributed  to  detecto  '  and  circuit  response  times. 

Since  the  detector  properties  were  only  approximately  known,  it 
was  difficult  to  determine  a  precise  value  for  the  peak  power  of  the 
idler.  The  best  estimate  (based  upon  the  detector  responsivity)  was 
that  the  observed  idler  power  was  approximately  5  W. 

While  the  experimental  results  are  important  taken  by  themselves, 
their  major  significance  for  this  work  are  in  their  comparison  with 
the  calculations  presented  in  the  previous  chapters.  Along  this  line, 
uhe  following  points  are  significant: 

(1)  The  experimental  results  confirm  that,  for  near- forward 
scattering,  the  gain  for  the  248  cm"1  mode  is  significantly  higher 
than  that  of  the  628  cm  ^  mode,  since  tunable  near- forward  scattering 
from  the  628  cm  1  mode  was  not  observed. 

(2)  For  the  case  of  normal  pimp  incidence,  the  experimentally 

inferred  idler  frequency  is  50  cm"1.  Referring  to  Fig.  4.5  it  is 

seen  that  the  calculations  show  g  to  have  a  maximum  at  v.  <=»  5°  cm  1 

s  1 

For  other  than  normal  incidence,  the  frequency  of  the  signal  output 
is  determined  primarily  by  the  feedback  provided  by  the  crystal  surfaces 

( 3)  The  calculated  value  of  g  for  Vj  *  50  cm  1  and  a  pump 

Q  p  S  ^  ^ 

intensity  of  5  x  10  W/cm  is  about  24  era"1  (see  Fig.  4.5).  If  the 
transverse  intensity  distribution  of  the  Gaussian  pump  beam  had  been 
included  in  the  theory,  the  calculated  gain  would  be  reduced  by  about 
a  factor  of  2.  In  addition,  the  effect  of  50$  pump  depletion  can 
be  included  in  a  very  approximate  manner  by  using  a  constant  pump 
intensity  of  three-quarters  of  its  input  value.  This  would  reduce 
the  calculated  gain  by  about  2S$.  Also,  surface  reflections  should 
be  taken  into  account.  The  net  effect  of  including  these  factors 

is  to  reduce  g  to  about  8  cm"  ,  which  is  in  reasonable  agreement 

s  -1 
with  the  experimentally  determined  value  of  5  cm 
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(4)  The  tuning  curve  agrees  with  that  predicted  by  the 
intersection  of  the  phase-matching  lines  with  the  theoretical 
dispersion  curve.  A  more  precise  confirmation  of  this  agreement 
is  provided  by  the  spontaneous  scattering  experiments  described 

later . 

(5)  The  temporal  response  pictures  point  out  two  significant 
facts.  First,  a  complete  theory  clearly  must  include  pump  depletion. 
Second,  the  fact  that  the  depletion  of  the  pump  occurs  within  the 
first  few  nanoseconds  of  the  pulse  indicates  that  the  use  of  a  steady 
state  theory  (for  the  remainder  of  the  pulse)  is  reasonable. 

(6)  The  detection  of  the  idler  is  of  itself  significant 
theoretically,  since  this  represented  the  first  confirmation  of  its 
existence  in  near-forward  scattering  from  vibrational  modes.  In 
addition,  the  measured  idler  power  of  ~5  W  is  consistent  with  the 
predicted  power  of  about  60  W,  (see  Fig.  4.6),  when  reflections  at 
the  crystal-air-interface  (~50$)  and  absorption  in  the  crystal  path 
length  from  the  generating  volume  to  the  crystal  surface  are  taken 
into  account. 

In  summary,  the  experimental  observations  were  found  to  be  in 
good  agreement  with  the  calculations. 

D.  SPONTANEOUS  SCATTERING  FROM  THE  248  cm"1  MODE. 

The  final  experiments  which  were  performed  involved  measurements 
of  the  optical  (signal)  power  which  is  generated  in  spontaneous 
scattering  from  the  248  cm  1  mode.  The  purpose  of  this  effort  was 
two-fold: 

(1)  Accurate  measurements  of  the  signal  frequency  as  a  function 
of  scattering  angle  could  be  made  to  verify  more  precisely  that  the 
tuning  does  follow  the  theoretical  dispersion  curve. 

(2)  The  variation  in  the  magnitude  of  the  scattered  signal 
power  as  a  function  of  idler  frequency  could  be  ascertained  and 
compared  to  calculations. 

The  experimental  set-up  was  the  same  as  that  described 
previously  (see  Fig.  5.4),  with  the  following  exceptions: 

a)  The  crystal  was  3  cm  long  with  polished  front  and  back 
faces . 
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b)  Standard  optical  gratings  were  used  in  the  Jarrel-Ash 
monochromators,  and  two  of  the  latter  were  used  in  tandem.  This 
reduced  the  monochromator  bandwidth  to  ~5  X  and  allowed  for  the 
detection  of  weak  signals  at  frequencies  close  to  a  strong  pump. 

(3)  The  detector  was  an  RCA  7200  photomultiplier  tube. 

The  results  of  these  experiments  are  shown  in  Figs.  5-12, 

5.13*  and  5.14.  In  Fig.  5*12,  measured  points  are  plotted  on  the 
dispersion  curve  described  in  Chapter  IV.  Within  experimental  error, 
the  measured  points  fall  directly  on  the  calculated  curve. 

The  agreement  with  theory  in  this  case  is  better  than  that 
obtained  in  the  stimulated  experiment.  This  lb  simply  because 
the  frequency  measurements  are  more  accurate  in  the  spontaneous 
experiments,  due  to  the  more  precise  measuring  techniques  which 
were  used. 

The  smallest  measureable  pump-to-s ignal  frequency  shift  was 
about  kO  cm  1  (AX  **  10  X),  corresponding  to  6  =  0.22°  .  This  was 

limited  by  the  presence  of  the  strong  pump  signal. 

The  output  of  the  detection  system  was  plotted  on  a  chart 
recorder  for  each  setting  of  the  angle  0  .  These  are  shown  in 

Fig.  5*13,  where  all  of  the  chart  records  have  been  drawn  in  one 
sketch.  In  this  case,  meaningful  data  could  be  obtained  only  for 
frequency  shifts  greater  than  about  110  cm-1  ( 25  X) .* 

The  magnitude  of  the  scattered  signal  (taken  from  the  chart 
recordings)  is  plotted  in  Fig.  5.1k,  along  with  the  calculated 
curve  (redrawn  from  Fig.  k. 7).  The  data  are  normalized  so  that 
the  measured  point  A  lies  on  the  calculated  curve.  The  agreement 
between  measurement  and  calculation  is  again  quite  good. 

In  summary,  therefore,  it  can  be  said  that  the  theory  for 
scattering  from  optical  modes  developed  in  Chapter  III  and  the 
experimental  results  described  in  Chapter  V  are  in  good  agreement. 

This  indicates  that  the  process  of  tunable  signal  and  idler  generation 
in  scattering  from  crystal  optical  modes  is  reasonably  described  by 
the  material  in  this  report. 


to  £0 

With  a  recently  reported  detection  technique,  it  may  be  possible 
to  extend  the  data  to  within  25  cm“^  of  a  strong  p\mp. 
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FIG.  5.1**--  Magnitude  of  scattered  aignal  power  as  a  function  of  idler 
frequency.  The  dots  are  the  measured  points  from  Fig.  5.15;  the  curve  is 
a  reproduction  of  a  portion  of  the  calculated  curve  of  Fig.  I*. 7.  The  data 
are  normalized  so  that  the  measured  point  "A"  lies  on  the  calculated  curve. 


CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  contributions  of  this  work  are: 

(1)  Development  of  the  theory  for  tunable,  near- forward, 
stimulated  and  spontaneous  light  scattering  and  the  application  of 
the  theory  to  the  transverse  modes  in  LiNbO^. 

(2)  Observation  of  tunable,  stimulated  optical  and  infrared 
outputs  in  scattering  from  the  248  cm"  mode  in  LiNbO^,  with  experi¬ 
mental  parameters  in  close  agreement  with  predicted  values. 

(3)  Construction  of  an  experimental  test  set  for  performing 
spontaneous  measurements,  consisting  of  an  argon  laser  and  a  phase- 
sensitive  detection  apparatus. 

(4)  Measurements  of  spontaneous  scattering  from  the  248  cm  ^ 
mode,  accurately  confirming  both  the  predicted  tuning  characteristic 
in  accordance  with  the  infrared  dispersion  curve,  and  the  calculation 
of  the  frequency  dependence  of  the  magnitude  of  the  scattered  signal. 

With  reference  to  further  extensions  of  this  work,  the  following 
suggestions  are  offered.  The  inclusion  of  the  effects  of  pump  deple¬ 
tion  and  transverse  intensity  variations  of  the  beams  would  improve 
the  validity  of  the  calculations.  On  the  experimental  side,  the 
characteristics  of  the  idler  radiation  generated  in  stimulated 
scattering  from  the  248  cm  1  mode  should  be  examined  in  more  detail. 
In  particular,  better  measurements  of  the  power  output  should  be 
made,  along  with  studies  of  the  linevidth  of  the  radiation.  In 
addition,  the  evaluation  of  other  crystals  will,  without  doubt, 
continued  to  play  a  central  role  in  the  field  of  nonlinear  optics. 
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APPENDIX  A 


THE  ENERGY  DENSITY  FUNCTION 


In  the  energy  density  approach,  the  existence  of  a  phenomeno¬ 
logical  energy  density  function  (U)  is  assumed,  and  a  method  for 
directly  deriving  the  driving  terms  (polarization  and  forces)  is 
given.  In  particular,  U  is  assumed  to  be  a  function  of  the  lattice 
displacement  coordinates  Q  and  the  electric  fields  E  ;  that  is, 

U  =  U(Q,E) 


and  the  function  is  defined  such  that  the  polarization  at  frequency 
ai.  ,  P(uu)  ,  and  the  force  at  frequency  ui^  ,  F(e^)  ,  are  given 

by: 


dU 

P(u>.) 

“  * 

1 

dE(a>±) 

1  dU 

F(o)  ) 

- - 

1 

N  SQ(u>. 

(A.l) 

y 


where  N  is  the  number  of  unit  cells  per  unit  volume.  The  dimen¬ 
sions  of  U  are  energy  per  unit  volume. 

To  illustrate  the  approach,  we  examine  the  linear  problem, 
where  only  a  vibration  ind  an  electric  field  at  u>^  are  considered. 

Since  the  polarization  und  the  force  will  be  linear  in  Q  and  E  , 

(2) 

we  begin  with  a  second  order  U  (written  as  U  3,  which  is  quadratic 
ir.  Q  and  E  .  Using  a  slightly  modified  form  of  Kleirman's  nota- 

25 

tion,  we  get: 

U(2)(Q,K)  -  -  ~  (Nb^Q2*  2Nfc12<|B  ♦ 


-  9e  - 


Writing  Q  and 


E  as 


Q  =  Q(to. )  +  Q(w.)*  >  E  =  E(o>  )  +  E(cd.  )* 

1  -L  11 

and  retaining  only  do  terms  in  *  We 

=  -  [Nbj.jQftV  b22EE*+  Nbl2(QVE+E*Q)] 

and  using  (A.l), 

P(u>.)  =  Nbl2Q  ♦  b22E 
*<*<)  «  bnQ  +  b12E 


Note  that  this  yields  the  expected  result,  xhe  polarization 
is  given  by  two  terms-one  due  to  the  vibration  and  the  other  being 
the  frequency  independent  term.  Similarly,  the  force  is  given  by 
a  term  proportional  to  Q  (the  restoring  force)  and  a  term  propor¬ 
tional  to  the  electric  field  E 

(5) 

To  examine  the  first-order  nonlinearity,  we  consider  U  ' 
which  is  cubic  in  Q  and  E  .  This  is  given  by 

U5(Q,E)  =  -  ±  [f^*  3f2Q h  +  3f5QE‘+  fuE?j  .  (A. 2) 

For  the  case  -e  ere  considering  in  the  text,  there  are  electric  fields 
at  a>  ,  a  and  a>,  ,  and  a  lattice  displacement  at  co.  .  Therefore, 

f  S  i  i 

we  write 

Q( lattice)  ■  Q(uj.)  ♦  Q(u>. )  ,  E(*vanp)  »  E(u>  )  +  E{as  ) 
ll  p  p 

E(  idler)  -  E(o>.  )  ♦  E(tu. )  ,  E(  signal)  «  E(u>  )+E(u>  ) 
XI  s  s 

By  dc  terms,  it  is  mgant  terms  which  exhibit  no  frequency  de¬ 
pendence.  It  can  be  shevrr0  that  only  tnese  terns  in  U  lead  to 
driving  terms  which  can  produce  significant  interaction. 
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Examining  (A. 2),  we  note  that  the  de  terms  can  arise  only  from 
the  Off  and  the  E^  terms.  Therefore,  reduces  to 

=  -  f,[Q(m.)  E{cd  )*e(cjd  )+  c.c.] 
j  i  s  p 

.  *  * 

- [E(o\)  E(cog)  E(co  )+  c.c.] 

■5  1  s  P 

Defining  f ^  s  d^  and  f^/3  =  NdQ  ,  we  arrive  at 

V  J  =  -  [d _E(<d .  )E(oj  )  E(co  )  +  Nd  E(o>  )E(o>  )  Q(o>.)  ]  +  c.c. 

ilipb  1  P  S  1 

(A. 3) 

Equation  (3-5)  is  an  extension  of  (A. 3}  to  a  multimode  material 


APPENDIX  B 

CALCULATION  OF  PHASE  MATCHED  GAIN  IN  "HIGH  LOSS" 

REGION  FOR  MULTI-MODE  MEDIA 

2  2  2 

Setting  g.  =  0  and  k  =  (ai /c  )  e  (i.e.,  phase  matching 

1  S  <3  00  g 

I’or  the  signal),  we  find  that  Eq.  (J.l4)  becomes 


C  (B.l) 

In  the  region  where  ox  is  outside  the  linewidth  of  each 
of  the  d  modes, 

*  2  2  2  2 
D_.  (o3.  )  =  Oi.  -  Ci),  •>  i£D.r.  ft  0)rt  -  CD. 

d  i  o,  i  i  i  o  i 
3  d 

Also,  the  imaginary  part  of  the  second  term  on  the  RHS  of 
(B.l)  is  small  and  can  be  neglected.  Thus,  we  must  consider  only 
the  imaginary  part  of  the  first  term.  This,  in  turn,  is  due  to  the 
imaginary  part  of  its  denominator.  Expanding  the  denominator  gives 


(B.2) 
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where  the  condition  that  (o>^,k^)  lie  on  the  undamped  dispersion 
curve,  that  is, 


has  been  imposed. 

Substitution  of  (B.2)  into  (B.l)  yields  the  final  result: 


APPENDIX  C 


SOLUTION  OF  BOUNDARY  VALUE  PROBLEM  TO 
DETERMINE  DIRECTION  OF  ~y 

The  direction  of  "7  can  be  determined  from  an  analysis 

of  the  boundary  conditions  which  must  be  satisfied  at  the 

interface  between  the  linear  (air)  and  nonlinear  media.  The 

4i 

approach,  which  was  suggested  by  Bloembergen  and  Pershan, 
makes  use  of  the  fact  that  a  necessary  and  sufficient  condition 
for  the  requirement  that  the  tangential  components  of  the  electric 
fields  be  continuous  across  a  boundary  is  that  the  components  of 
the  momentum  wave  vector  (It)  parallel  to  the  boundary  remain 
conserved. 

This  is  now  applied  to  the  geometry  of  Fig.  C.l.  The 

boundary  is  the  y  =  0  plane,  and  the  pump  wave  vector  external 

to  the  nonlinear  medium  (It6  )  makes  an  angle  0  with  the  normal 

P 

to  the  boundary.  The  phar  -matching  triangle  inside  the  material 
is  also  shown. 

For  the  purposes  of  this  analysis,  the  space-dependent  part 
of  plane  wave  expressions  for  the  signal  and  idler  will  be  written 

it-r  iO?'-i£")-? 
e  =  e  * 


x 


FIG.  C.l--Sketch  of  the  geometry  for  the  boundary  value  problem. 
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where  k  =  y  ,  used  in  the  text.  Continuity  of  the  x-component 
wave  vector  leads  to 


k‘x  =  k*x  =  k*x  +  k.  x 
P  P  SI 


(C.l) 


r»e 


where  x  is  a  unit  vector  in  the  x-direction.  Since  k  is 

P 

real,  (C.l)  yields 


x  = 


k  •  x  +  k.  *  x 
s  1 


k  •  x  =  -  k.  •  x 

s  1 


But,  it  has  been  shown  [Eq.(?.llb)]  that 

•  -7*/  /  A  t>/  /  A  A 

..  k  •  x  =  k  •  x  =  0  . 

s  .  ' 

1 

or  "7  lies  along  the  normal  to  the  boundary  (the  y  axis  in  the 
geometry  of  Fig.  C.l). 

Recalling  the  definitions  of  g.  and  g  in  Eq.  (3.13),  we 

1  S 

obtain  immediately  the  desired  relationship: 

«i  -  «.  •  <c-2> 

For  near  forward  scattering  and  for  near-normal  incidence 
of  the  pump,  Eq.  (C.2)  reduces  to: 


g.  =  gg  cos  cp 


(C .  3) 
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APPENDIX  D 

SOLUTION  FOR  g  IN  NONCOLLINEAR  LOW  LOSS  REGION 
s 


Setting  k2  =  (oo2/  c2)  e 


s  ~  's'  '  «  in  the  secular  equation  (3.1*0 


and  substitution  of  g  =  g./cos  <P  yields: 

S  1  .2 


l-k.  -  ik 

X 


®i  *\ 

a  +  ?  eiJ 


-§  (|AJ2)  E-Sf-^  +  ikA 


S  1 


2  2 
05.05 
1  S 


where 


p  v5  j  d3h >7 


(|Ap|2)  =  |Ap|2  cos  9 


(D.l) 


Thus,  g  has  been  eliminated  from  the  equation  and  proceeding 
s 


as  in  Appendix  B  of  HG  gives: 

a 

g  •  “  ~ 

1  2 


where 


1 6  a'2V/2 


1  + 


-  1 


a. 

i 


>2.  2 

a'  =  a  cos  m 
P  P  y 


(D.2) 


and  therefore 


cos  cp 


(D.3) 
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APPENDIX  E 

SCATTERED  IDLER  POWER 


The  idler  power  scattered  into  a  small  solid  angle  and 
frequency  range  is  given  by: 


Pi} 


a a. 

i  i 


lP  dQ.do). 
p  1  1/ 


P 

P  i 


(E.l) 


where  the  prime  indicates  that  the  solid  angle  is  measured  out¬ 
side  the  crystal  (after  refraction  at  the  crystal  urface) . 

For  an  aperture  of  height  H  located  at  distance  r  from 
the  end  face  of  the  crystal,  the  accepted  solid  angle  is 

(see  Fig.  E.l) : 


AB  »  -g 
r 


HrAcp 


HAcp 


and 


Pi> 


P. 

1 


aa.'do). 

p  i 


H*q> 

j  P  - An). 


(E.2) 


The  expressions  derived  in  the  text,  however,  give  the 
scattered  power  per  unit  solid  angle,  with  the  solid  angle  measured 
about  the  signal  propagation  direction.  To  convert  this  to  the  form 
required  for  Eq.  (E.l),  we  write: 


But 


(E.3) 
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FIG.  E.l — Sketch  of  geometry  for  the  solid  angle  calculation. 
Therefore, 


P^/nyv 


d£  Vr/Vi. 
de 


(E.4) 


which  is  the  desired  result. 

It  should  he  recalled  that  specifying  the  solid  angle 
specifies  the  frequency  range  through  the  intersections  of 

the  phase-matching  lines  with  the  dispersion  curve. 
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